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Abstract 

We discuss the unification of gauge couplings within the framework of 
a wide class of realistic free-fermionic string models which have appeared 
in the literature, including the flipped SU(5), 50(6) x 50(4), and various 
SU(3) x SU(2) x U(l) models. If the matter spectrum below the string scale 
is that of the Minimal Supersymmetric Standard Model (MSSM), then string 
unification is in disagreement with experiment. We therefore examine several 
effects that may modify the minimal string predictions. First, we develop a 
systematic procedure for evaluating the one-loop heavy string threshold cor- 
rections in free-fermionic string models, and we explicitly evaluate these cor- 
rections for each of the realistic models. We find that these string thresh- 
old corrections are small, and we provide general arguments explaining why 
such threshold corrections are suppressed in string theory. Thus heavy thresh- 
olds cannot resolve the disagreement with experiment. We also study the 
effect of non-standard hypercharge normalizations, light SUSY thresholds, and 
intermediate-scale gauge structure, and similarly conclude that these effects 
cannot resolve the disagreement with low-energy data. Finally, we examine 
the effects of additional color triplets and electroweak doublets beyond the 
MSSM. Although not required in ordinary grand unification scenarios, such 
states generically appear within the context of certain realistic free-fermionic 
string models. We show that if these states exist at the appropriate thresholds, 
then the gauge couplings will indeed unify at the string scale. Thus, within 
these string models, string unification can be in agreement with low-energy 
data. 



* E-mail address: dienes@sns.ias.edu 
t E-mail address: faraggi@sns.ias.edu 



1 Introduction 



LEP precision data provides remarkable confirmation of the Standard Model of 
particle physics. However, many fundamental problems are not addressed in the 
context of the Standard Model, leading to the expectation that a more fundamental 
theory must exist in which the Standard Model appears as an effective low-energy 
limit. While many possible extensions of the Standard Model are highly constrained 
or ruled out by experiment, supersymmetric theories are in agreement with all avail- 
able data. Furthermore, the top-quark mass range required in supersymmetric sce- 
narios of electroweak symmetry breaking is in agreement with that suggested by 
CDF/DO direct observation J2J and LEP precision data ||. In recent years it has 
also been suggested that the success of gauge coupling unification in the Minimal 
Supersymmetric Standard Model (MSSM) [[| provides evidence for the validity of 
supersymmetric grand unified theories (SUSY GUT's). 

While SUSY GUT's provide a useful parametrization of the sparticle spectrum at 
low energies and of the boundary conditions at the GUT scale, they are incomplete 
theories. First, they do not explain the origin of the Standard Model spectrum. 
Second, in order to evade proton-lifetime constraints, some ad hoc global symmetries 
must be imposed, and a doublet-triplet splitting mechanism is required. Finally, 
despite the fact that the unification scale is just one or two orders of magnitude 
below the Planck scale, a consistent treatment of quantum gravity is lacking. 

Remarkably, superstring theories [EJ , the only available candidates for a consistent 
theory of quantum gravity, accommodate N = 1 supersymmetric theories as their 
low-energy effective theories. Even more remarkably, heterotic string theories || 
provide a general framework in which the origin of the observed particle spectrum 
and interactions may be understood |7] and the doublet-triplet splitting problem may 
be resolved [§]. 

Among the string models constructed to date, the most phenomenologically realis- 
tic have been formulated within the so-called "free-fermionic" construction [^|, O]. 
Indeed, within this construction, many three-generation models can be obtained. 
While this may be an accident, it is more likely to be a reflection of some fundamental 
properties of string compactification. The free-fermionic construction is formulated 
at a highly symmetric point in the compactification space at which spacetime sym- 
metries are maximally enhanced. It turns out that the realistic free-fermionic models 
admit a Z2 x Z2 orbifold structure with standard embedding which is realized in this 



construction through the so-called "NAHE set" |Tj] of fermionic boundary-condition 
basis vectors. Such Z2 x Z2 orbifolds possess a structure which can naturally accom- 
modate three generations due to the existence of exactly three twisted sectors. These 
result from the action of the Z 2 x Z 2 twist on a six-dimensional compactified space. 
In general, the Z2 x Z2 orbifold does not produce a number of fixed points that can 
be reduced to three generations. However, precisely at the free-fermionic point in the 
toroidal compactification space, the number of fixed points is such that a reduction 
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to three generations can be achieved. Indeed, the number of fixed points from each 
twisted sector is reduced to one. Thus, each twisted sector produces one of the chiral 
generations of the Standard Model. 

The realistic free-fermionic models achieve remarkable success in trying to explain 
the different features of the Standard Model spectrum. In addition to naturally 
producing three generations with Standard Model gauge group, they also provide a 
plausible explanation for the heavy top-quark mass and the fermion mass hierarchy. 
Indeed, the realistic standard-like models suggest that only the top quark mass term is 
obtained at the cubic level of the superpotential, while the lighter fermion mass terms 
are obtained from non-renormalizable terms which are naturally suppressed relative 
to the leading cubic- level mass term. In fact, in this way a successful prediction of the 
mass of the top quark was obtained in Ref . [13] , three years prior to its experimental 
observation. Furthermore, an analysis of non-renormalizable terms up to eighth order 



then reveals how the fermion masses and mixing angles may also be generated (14 ] . 

The realistic free-fermionic models are therefore very appealing from a theoretical 
point of view, and may successfully explain the different features of the Standard 
Model. However, string theory in general, and the free-fermionic models in particu- 
lar, predict that the string unification scale is related to the Planck scale, and should 
be numerically of the order of (9(5 -10 17 ) GeV. Thus, a factor of approximately twenty 
separates the string unification scale from the usual MSSM unification scale extrap- 
olated from low-energy data. This discrepancy is one of major problems confronting 
string model-building. 

One possible solution is to construct string GUT models in which the GUT sym- 
metry is broken in the effective low-energy field theory [|1|, However, no realistic 
models of this sort have yet been constructed. Moreover, in this scenario, the prob- 
lems with proton-lifetime constraints reemerge. Indeed, this problem is more severe in 
string GUT's due to the possible appearance of baryon- and lepton-number- violating 
dimension-four operators |J and the anticipated difficulty in implementing the GUT 
doublet-triple splitting mechanism in string models. 

Another solution is provided if additional thresholds exist in the desert between 
the electroweak scale and the string unification scale. In fact, the availability of 
such additional thresholds in realistic free-fermionic models has been demonstrated 



in Ref. [|17]]. In this respect, imposing the restriction that the spectrum below the 
string scale is just that of the MSSM is ad hoc, and may be too restrictive. The 
successful unification of gauge couplings within the MSSM would then indicate that 
the MSSM is only an approximation to the complete theory between the weak scale 
and the Planck scale. 



A third suggestion, due to Ibanez fll8fl ) is that in string models the normalization 
of the weak hypercharge is, in general, different from that in grand-unified theories, 
and may just have the right value to allow string unification, even if the spectrum 
below the string scale is that of the MSSM. It is found that if k\ is in the range 
1.2 < k\ < 1.4, then unification at the string scale can be consistent with low-energy 
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data. However, whether such values of k\ can be achieved within consistent string 
models remains an open question. Moreover, we will argue that in string models we 
must have k\ > 5/3. 

A final possibility is that the contributions arising from the infinite towers of 
heavy string states can modify the above string tree-level predictions. Within this 
scenario, one would hope that these heavy string states can give rise to additional 
heavy string threshold corrections which shift the unification point down to (9(10 16 ) 
GeV. 

In this paper, we systematically re-examine the problem of gauge coupling uni- 
fication within the context of a wide variety of realistic free-fermionic models which 
have appeared in the literature. There are various reasons why this is an important 
undertaking. For example, while many of the possible effects that we consider have 
been discussed previously, each was treated in an abstract setting and in isolation. 
However, within the tight constraints of a given realistic string model, the mecha- 
nisms giving rise to three generations and the MSSM gauge group may ultimately 
prove inconsistent with, for example, large threshold corrections or extra non-MSSM 
matter. Moreover, the increased complexity of the known realistic string models may 
substantially alter previous expectations based on simplified or idealized scenarios. 
It is therefore important to rigorously calculate all of these effects simultaneously, 
within the context of a wide variety of actual realistic string models, in order to 
determine which path to successful gauge coupling unification (if any) such models 
actually take. 

Our analysis proceeds in several steps. First, we study the role of threshold correc- 
tions due to the infinite tower of heavy string states. In this we follow the definition 
of threshold corrections as given by Kaplunovsky [|19| , and explicitly calculate these 
corrections for a wide range of realistic free-fermionic models including the flipped 
577(5) "revamped" model ||, several SU(3) x SU(2) x U(l) models gl], ||, g|], 
an SO (6) x SO (4) model and even models with non-supersymmetric spacetime 
spectra. Our evaluation of these threshold corrections is done in two stages. First, 
we analytically evaluate the traces over the entire Fock space corresponding to a 
given string model, with various relevant combinations of gauge charge operators in- 
serted into the trace. Then, for each combination of gauge charge insertions, these 
results are expanded level-by-level, and the contributions from the states at each 
string energy level are integrated over the modular-group fundamental domain. In 
this manner, we can obtain results to any desired accuracy. Moreover, as we shall 
see, various non-trivial consistency checks can be performed. 

Our results show that threshold corrections due to the massive string states are 
small in free-fermionic string models. This result is a priori surprising, given the 
infinite numbers of heavy string states which potentially contribute to the string 
threshold corrections, but we are able to provide a general argument which explains 
why the threshold corrections from the massive string states are naturally suppressed 
in string theory. Our argument also explains why threshold corrections can grow large 
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only for large values of the string moduli. 

Given these results, we then proceed to systematically examine the other ef- 
fects which might potentially alleviate the discrepancy between the GUT and string 
scales. As discussed above, these include the effects of stringy non-standard U(l) hy- 
percharge normalizations, light SUSY thresholds, intermediate gauge structure, and 
additional matter beyond that predicted by the MSSM. We find that the effects of 
hypercharge normalizations, light SUSY thresholds, and intermediate gauge struc- 
ture are not sufficient to remove the discrepancy. By contrast, we surprisingly find 
that only the presence of additional matter in these models (in particular, certain 
color triplets and electroweak doublets with special hypercharge assignments) has 
a profound effect on the running of the gauge couplings, and precisely this matter 
appears naturally in a variety of the realistic string models. Indeed, within these 
models, we show the gauge couplings can indeed unify at the string scale when all 
of the above effects are taken into account. Thus, for such models, the disagreement 
between the GUT scale and the string scale can be naturally resolved. 

It is remarkable that string theory, which predicts an unexpectedly high unifi- 
cation scale M string , in many cases compensates by simultaneously also predicting 
precisely the extra exotic particles needed to reconcile this higher scale with low- 
energy data. Moreover, our analysis shows that while some string models naturally 
contain the extra matter needed to resolve this disagreement, other models do not 
and can actually be ruled out on this basis. Thus, the appearance of such extra 
matter becomes a low-energy prediction of these models which may be accessible to 
present-day experiments. 

We stress that this is the first time that such an exhaustive examination of all 
possible effects has been performed within the context of actual realistic string mod- 
els, and within the constraints that these models impose. Hence, one of the main 
results of our analysis is the new observation that only the appearance of extra exotic 
matter in particular representations can possibly resolve the experimental discrepan- 
cies. Indeed, because we have been able to rule out all other possible within these 
models, the appearance of such extra exotic matter becomes a prediction of successful 
string-scale unification. It is of course an old idea that the presence of extra matter 
can resolve the discrepancy between the GUT and string unification scales. What 
is highly non-trivial, however, is that this now appears to be the only way in which 
realistic string models can solve the problem. 

This paper is organized as follows. In Sect. 2 we review the general features of the 
realistic free-fermionic models, and in Sect. 3 we discuss how heavy string threshold 
corrections may be evaluated within the context of such models. In Sect. 4 we then 
explicitly calculate the threshold corrections within the flipped £77(5), £0(6) x £0(4), 
and SU(3) x SU(2) x £7(1) string models, and in Sect. 5 we test whether these 
models are in agreement with low-energy data by systematically analyzing these 
results along with the effects due to light SUSY thresholds, intermediate gauge and 
matter thresholds, as well as two-loop and Yukawa-coupling effects. It is here that 
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we find, contrary to naive expectations based on comparing Mqut and M str ing, that 
certain string models may be in full agreement with low-energy data and grand 
unification. Sect. 6 then contains a general argument based on the modular properties 
of the threshold corrections which explains why they must be small in the free- 
fermionic models, and Sect. 7 contains our conclusions. Finally, an Appendix contains 
an explicit listing of all of the string models that we will be considering in this paper. 
A short summary of some of the main results of this paper can also be found in 



Ref. fl25f 



2 Realistic free-fermionic models 

In the free-fermionic formulation of the heterotic string ||, all of the worldsheet 
degrees of freedom needed to cancel the conformal anomaly are represented in terms 
of internal free fermions propagating on the string worldsheet. In four dimensions, 
this requires 44 real left-moving (Majorana-Weyl) fermions and 20 real right-moving 
Majorana-Weyl fermions (or equivalently, half as many complex fermions, or any 
such consistent combinations of real and complex fermions). Under parallel trans- 
port around a non-contractible loop on the toroidal string worldsheet, these fermionic 
fields can generally accrue a phase, and each set of specified phases for all worldsheet 
fermions around all such non-contractible loops is called the "spin structure" of the 
model. Such spin structures are usually given in the form of boundary-condition 
"vectors", with each element of the vector specifying the phase of a corresponding 
worldsheet fermion. The possible spin structures which can be used in the construc- 
tion of string models are constrained by various string consistency requirements (e.g., 
the existence of a proper worldsheet supercurrent, proper spacetime spin-statistics 
assignments, physically sensible projections, and modular invariance). A model is 
therefore constructed by choosing a set of boundary condition vectors which satisfy 
these constraints. In general, these basis vectors must span a finite additive group 
H = J2k n k)°k where = 0, • • • , N Zk — 1. The physical massless states in the Hilbert 
space of a given sector a G H are then obtained by acting on the vacuum state of 
that sector with the worldsheet bosonic and fermionic mode operators, and by sub- 
sequently applying the generalized GSO projections. The U(l) charges Q(/) with 
respect to the unbroken Cartan generators of the four-dimensional gauge group are in 
one-to-one correspondence with the U(l) currents /*/ for each complex worldsheet 
fermion /, and are given by: 

QUO = \<*U) + F(f). (2.i) 

Here a(f) is the boundary condition of the worldsheet fermion / in the sector a; 
at(f) = for Neveu-Schwarz boundary conditions, and oc(f) = 1 for Ramond. Like- 
wise, F a (f) is a fermion-number operator counting +1 for each mode of / (and —1 for 
each mode of /*, if / is complex). For periodic complex fermions [i.e., for a(f) = 1], 



6 



the vacuum is a spinor with two degenerate vacuum states |+) and |— ). These states 
are respectively annihilated by the zero modes /o and /o* which obey the Clifford 
algebra, and have fermion numbers F(f) = 0, —1 respectively. 

In the realistic free-fermion models that we will be considering, the worldsheet 
fermions are as follows: 

• a complex right-moving fermion, denoted ip 11 , formed from the two real fermionic 
superpartners of the coordinate boson X M ; 

• six real right-moving fermions denoted x '"' 6 > often paired to form three com- 
plex right-moving fermions denoted x 12 > X 34 > an d X 56 > 

• 12 real right-moving fermions, denoted y 1 '-' 6 and u; 1 ' "' 6 ; 

• 12 corresponding real left-moving fermions, denoted y 1 '-'® and 57 1 '-' 6 ; and 

— 1 5 

• 16 remaining complex left-moving fermions, denoted ip ' "' , rj 



T ' ,3 , and 



The realistic models in the free-fermionic formulation are then generated by spec- 

, WH for 



ifying a special basis of boundary-condition vectors |20|, [26|, |24], |21], [13], ^2| 
these worldsheet fermions. This basis is constructed in two stages. The first stage 
consists of introducing the so-called NAHE set f22|, which is a set of five boundary 
condition basis vectors denoted {1, S, hi, b 2 , b 3 }. With '0' indicating Neveu-Schwarz 
boundary conditions and '1' indicating Ramond boundary conditions, these vectors 
are as follows: 
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(2.2) 



As can be seen, the vector 1 has periodic boundary conditions for all the worldsheet 
fermions. The vector S generates the spacetime supersymmetry, and the sectors bi, 
b 2 , and b 3 correspond to the three twisted sectors of a Z 2 x Z 2 orbifold. Corresponding 
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to this set of boundary-condition vectors is the following choice of phases which define 
how the generalized GSO projections are to be performed in each sector of the theory: 
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1 



(2.3) 



The remaining projection phases can be determined from those above through the 
self-consistency constraints. The precise rules governing the choices of such vectors 
and phases, as well as the procedures for generating the corresponding spacetime 
particle spectrum, are given in Refs. |], |10| . The mapping between the notation used 
here ]10[ and the notation used in Ref. || is given in the Appendix. 

After imposing the NAHE set, the resulting model has gauge group 50(10) x 
50 (6) 3 x Eg and N — 1 spacetime supersymmetry. The vector 5* is the supersymme- 
try generator and the superpartners of the states from any given sector a are obtained 
from the sector 5 + a. The spacetime vector bosons that generate the gauge group 
arise from the Neveu-Schwarz sector and from the sector / = 1 + bi + b 2 + b 3 . The 
Neveu-Schwarz sector produces the generators of 50(10) x 50(6) 3 x 5*0(16). The 
sector 1 + b x + b 2 + b 3 produces the spinorial 128 of 50(16) and completes the hid- 
den gauge group to Eg. The vectors bi, b 2 , and b 3 correspond to the three twisted 
sectors in the corresponding orbifold formulation and produce 48 spinorial 16's of 
5*0 (10), sixteen from each of the sectors b 1; b 2 , and b 3 . 

As can be seen from ( j2.2| ), the NAHE set divides the 44 left-moving and 20 
right-moving real internal fermions in the following way: the complex fermions 
■0 1 ''"' 5 produce the observable 5*0(10) symmetry; the complex fermions (j) 1 ''"' 8 pro- 



duce the hidden E 8 gauge group; and the fermions {77 1 , y 3, '"' 6 }, {r] 2 ,y 



1,2 



CO'' 



and 



{rj 3 , u; 1 ''"' 4 } give rise to the three horizontal 50(6) symmetries 
fermions {y, uj} are divided as well, in groups {y 3 ' "' 6 } 



The left-moving 



{y 1 ' 2 ,^ 5 ' 6 }, and K<-< 4 }. 
The left-moving fermions \ 12 i X M i X 56 carry the supersymmetry charges. Each sector 
bi, b 2 , and b 3 imposes periodic boundary conditions for the fermions {%p^\%p l, '" )5 ) and 
for (X12, {y 3 '"'' 6 \y 3 '"' 6 }, rj 1 ), (X34, {y 1 ' 2 , ^ 5 ' 6 |y 1 ' 2 ^ 5 ' 6 }, r) 2 ), or (xs6, {w 1 '-' 4 !^ 1 -- 4 }, r/ 3 ) re- 
spectively. This division of the internal fermions is a reflection of the equivalent 
underlying Z 2 x Z 2 orbifold compactification | 12| . The set of internal fermions 
{y, u\y, a)} 1 ''"' 6 corresponds to the left/right symmetric conformal field theory of the 
heterotic string, or equivalently to the six-dimensional compactified manifold in a 
bosonic formulation. This set of left / right symmetric internal fermions plays a fun- 
damental role in the determination of the low-energy properties of the realistic free- 
fermionic models. 

The second stage in the construction of the realistic models consists of adding 
three additional basis vectors to the above NAHE set. These three additional basis 
vectors, which are often called {«,/?, 7}, correspond to "Wilson lines" in the orbifold 
construction. The allowed fermion boundary conditions in these additional basis vec- 
tors are of course also constrained by the string consistency constraints, and must 
preserve modular invariance and worldsheet supersymmetry. The choice of these ad- 
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ditional basis vectors {a, (3, 7} nevertheless distinguishes between different models 
and determine their low-energy properties. For example, three additional vectors 
are needed to reduce the number of massless generations to three, one from each 
sector bi, b2, and b3, and the choice of their boundary conditions for the internal 
fermions {y, u\y, a)} 1 ''"' 6 also determines the Higgs doublet-triplet splitting and the 
Yukawa couplings. These low-energy phenomenological requirements therefore im- 
pose strong constraints p2| on the possible assignment of boundary conditions to the 
set of internal world-sheet fermions {y, u\y, a)} 1 ''"' 6 . 

One then finds a variety of possibilities, with the 50(10) gauge symmetry broken 
to one of its subgroups, either SU(5) x 17(1), 50(6) x 50(4) or SU(3) x SU{2) x 
U(1)b-l x U(1)t 3r - These breakings are achieved by the assignment of the following 
boundary conditions to the fermion set ip\/2- To achieve a breaking to SU(5) xU(l), 
we need to assign boundary conditions 1/2 to all of these fermions simultaneously. By 
contrast, to achieve a breaking to 50(6) x 50(4), we assign boundary condition 1 to 
only the first three fermions, while assigning boundary condition to the remaining 
two. Finally, to break the 50(10) symmetry to SU(3) x SU(2) x U(l) c x U(l) L , we 
impose both of these breakings via two separate basis vectors.^ The complete listings 
of boundary conditions and GSO phases for each of the models we will be considering 
in this paper can be found in the Appendix. 

Note that all of these realistic free-fermionic models have three Z7(l) symmetries, 
denoted by U(l) r , (j = 1,2,3); these are respectively generated by the left-moving 
world-sheet currents fjl/ 2 Vi*/2^ ^1/2^1/2' anc ^ ^1/2^1/2 • These three U(l) symmetries 
arise due to the Z2 x Z2 twist with standard embedding. Additional horizontal Z7(l) 
symmetries, denoted by U(l) rj (j = 4, 5, ...), arise by pairing two real fermions from 
the sets {y 3, "' ,e }, {y 1 ' 2 , a) 5 ' 6 }, and {a; 1 ''"' 4 }. The final observable gauge group depends 
on the number of such pairings. Indeed, the rank of the entire gauge group can vary 
between r = 16 to r = 22. To every one of the horizontal gauged right-moving £7(1) 
symmetries corresponds a horizontal left-moving global £7(1) symmetry. Finally, the 
realistic free-fermionic models also contain Ising-model operators that are obtained 



by pairing a right-moving real fermion with a left-moving real fermion [28 . 

The hidden sector in the free-fermionic standard-like models is determined by the 
boundary condition of the internal left-moving fermions 1 ' '"' 8 . In the NAHE set, 
the contribution to the hidden Eg gauge group comes from the Neveu-Schwarz sector 
and from the sector I = 1 + hi + b 2 + b 3 , which produces the adjoint and spinorial 
representations of 50(16) respectively. The final hidden gauge group is determined 
by the assignment of boundary conditions for the worldsheet fermions cf) 1 ' '"' 8 from the 
vectors {a,/?, 7}, and by the choices of generalized GSO projection phases. 

The realistic free-fermionic models contain three chiral generations from the sec- 
tors bi, b2, and b3. At the level of the NAHE set alone, there are 48 generations, 
16 from each of these sectors. However, the vector 7 reduces the number of gener- 



" Recall that U(l)c = %U(1) B -l and that U(1) L = 2U(1)t 3r . 
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ations by a factor of two by fixing the charge under U(1)r. for each of the sectors 
bi, b 2 , and b 3 to be either +1/2 or —1/2. A further reduction to three generations 
is then obtained by carefully choosing the {a, (3, 7}-boundary conditions for the real 
fermions {y, u\y, a)} 1 '"'' 6 . Each one of the vectors {a,/3, 7} reduces by a factor of 
two the number of degenerate vacua that arise from the real fermions {y, u\y, a)} 1 ''"' 6 
in each of the sectors bi, b2, and b3. After the GSO projections from the vectors 
{a, (3, 7}, each sector b 1; b 2 , and b 3 produces one chiral generation. Thus, the as- 
signment of boundary conditions for the set {y, u\y, cj} 1 ''"' 6 from the vectors {a, /3, 7} 
is constrained by requiring three light generations. Further constraints on the assign- 
ment of boundary conditions for the set {y, u\y, a)} 1, "' 6 from the vectors {a, j3, 7} are 
imposed by requiring other phenomenological criteria, such as the presence of Higgs 
doublets in the massless spectrum, the projection of colored Higgs triplets, the exis- 
tence of a phenomenologically realistic supersymmetric vacuum, and the existence of 
non- vanishing Yukawa couplings that may produce a realistic fermion mass spectrum. 
Satisfying all of these phenomenological criteria simultaneously is a highly non-trivial 
task, and it is indeed a remarkable feat that models which successfully incorporate 
all of these features have been constructed. 



3 Calculating Threshold Corrections in Free-Fermionic 
Models 

In this section we first provide a general model-independent summary of how 



threshold corrections are defined in string theory Hl9| . We then outline our specific 
procedure for calculating threshold corrections within the context of free-fermionic 
string models. 

Preliminary Comments Concerning String Unification and Mass Scales 

We begin, however, with some preliminary comments concerning the unification 
of gauge couplings and mass scales within string theory. 

In order to relate string dynamics to low-energy experimental observables, we are 
interested in the effective field theories of the light (massless) string excitations. In 
particular, we are interested in the effective gauge couplings of the four- dimensional 
gauge group, and their beta functions. Unlike the case in field theory, however, in 
string theory all couplings and mass scales are dynamical variables, and are ultimately 
set by the values of certain moduli fields. Consequently, determination of the gauge 
couplings, and especially their "unification scale" , becomes a highly non-trivial issue. 

Perhaps the most important factor in this regard is the presence of moduli fields 
in the string spectrum. Such moduli fields, which are massless gauge-neutral Lorentz 
scalar fields whose effective potential is classically and perturbatively flat, arise as 
follows. Recall that the consistency of a heterotic string theory requires 22 internal 
bosonic degrees of freedom in the non-supersymmetric sector (or more generally, an 
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internal conformal field theory of total central charge 22). Sixteen of these bosonic 
degrees of freedom are compactified on a flat torus with fixed radius, and generically 
give rise to the internal symmetries of the theory which are interpreted as the gauge 
symmetries of the low-energy field theory. The remaining six bosonic degrees of 
freedom in the non-supersymmetric sector, combined with six compactified bosonic 
degrees of freedom in the supersymmetric sector, are compactified on an internal 
manifold that can be a Calabi-Yau manifold or an orbifold. The moduli fields are 
then precisely those massless fields which parametrize the size and shape of this 
six-dimensional internal manifold. 

In general, the (left-moving) gauge symmetry that one obtains for generic points 
in the moduli space of four-dimensional strings is U(l) 22 . At special points in the 
moduli space, however, the gauge symmetries may be enhanced, with additional 
vector states becoming massless. Indeed, such points with enhanced symmetries 
can often be realized as the results of compactifications of higher- dimensional string 
theories, with the enhancement of the gauge symmetry arising for special choices of 
the compactification moduli (or equivalently, special choices of background fields) 
p9|] . In this paper, we are interested in such points at which the effective gauge 
symmetry is enlarged to contain at least the MSSM gauge group SU(3)c x SU{2) L x 
U(1)y- Similarly, as one moves around in moduli space, there also exist special 
points at which the internal (c r i g ht, Cieft) = (9, 22) worldsheet conformal field theory 
effectively decomposes into a product of smaller worldsheet conformal field theories. 
For example, in the class of models that we study in this paper, all of the internal 
degrees of freedom can be represented in terms of free fermions propagating on the 
string worldsheet. Hence, each of the phenomeno logical models that we study in this 
paper corresponds to a special set of fixed values for the moduli fields. 

The gauge couplings of the four- dimensional gauge group factors are all gener- 
ally related to the gauge coupling of the string theory in ten dimensions, and are 
consequently related to each other at tree level. Thus, a natural unification of gauge 
couplings occurs in string theory regardless of the presence of any grand-unified gauge 
group. At tree level, the size of these couplings at unification is determined |3(J by 
the VEV of a special moduli field, the dilaton (ft, via gi ~ e~^> . Unfortunately, the 
presence of a classically flat dilaton potential implies that one does not know, a pri- 
ori, the value of the dilaton VEV. Hence, at the classical level, one does not know 
the general size of the string coupling constant at unification. 

This observation is not just limited to the gauge couplings, but applies to the 
gravitational coupling as well. Indeed, because string theories naturally incorporate 
quantum gravity and contain a massless graviton in their spectra, the gauge cou- 
pling in the effective four- dimensional gauge group can ultimately be related to the 
gravitational coupling by calculating the scattering amplitude between two gauge 
bosons and the graviton |HJ. Since the string coupling constant is universal at the 



classical level (modulo Kac- Moody levels), and since this coupling is determined by 
the vacuum expectation value of the four-dimensional dilaton field, the scale for the 
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gravitational coupling is therefore also set by the value of the dilaton VEV. Thus all 
of the string couplings are related at tree level to a single string coupling constant 
whose value is a priori unknown. 

Even more importantly, the string "unification scale" M str i ng is also not a fixed 
quantity, for the specific choice of M str i n g is ultimately a matter of definition which 
depends on the chosen renormalization scheme. This will be discussed below. In 
general, this scale is related to the string tension (or the phenomenological Planck 
scale) via M string ~ \j\fa' ~ fi'string^pianck- It is an important feature that M str i ng 
and Mpi anc k be of roughly the same order of magnitude; otherwise, the effective four- 
dimensional theory will not be weakly coupled ]52" |, and our subsequent perturbative 
one-loop analysis will not be valid. 

The above considerations apply at tree level in string perturbation theory. How- 
ever, at the one-loop level, this situation becomes even more non-trivial, for in general 
the quantum corrections are sensitive to not only the massless string modes, but also 
the entire infinite tower of massive string modes, and these in turn generally depend 
critically on the expectation values of all of the moduli fields. In particular, the quan- 
tum corrections to the gauge couplings — i.e., the so-called "threshold corrections" 
to be discussed below — are complicated functions of the moduli fields, and hence the 
running of the gauge couplings of the four- dimensional group factors in the effective 
low-energy theory generally depends on the expectation values of all of these moduli 
fields. This important dependence has been studied by several groups |JJ| [53, [5^] , and 
is crucial for understanding how an effective potential might be generated which not 
only selects a preferred string vacuum (thereby lifting the degeneracy of string vacua 
corresponding to different values of the moduli VEV's), but which might also pro- 
vide a means of dynamical supersymmetry breaking (e.g., through the formation of 
gaugino condensates). Furthermore, understanding the dependence of the threshold 
corrections on the moduli fields gives great insight into the structure and symmetries 



of the string effective low-energy supergravity Lagrangian |34 . 

It should be noted, however, that most of these studies of the moduli dependence 
of the threshold corrections make some simplifying assumptions that do not hold for 
the realistic string models that we will be studying here. For example, the moduli de- 
pendence in those studies is usually extracted for the case of heterotic string theories 
with (2, 2) worldsheet supersymmetry^ and for restricted classes of toroidal compact- 
ifications. However, the models that we will be studying here correspond to (2, 0) 
compactifications in which the moduli fields are not simple to identify. Similarly, this 
class of models contains twisted moduli and sectors that correspond to Wilson lines 
in an orbifold formulation. Therefore, the previous results on the moduli dependence 
of the gauge couplings are in general not applicable to the realistic string models 
we will be studying, and we will need to develop a different strategy. Of course, as 
discussed above, in our setting the moduli fields are fixed at particular values in the 



* Recent results concerning the moduli-dependence of threshold corrections in (2,0) compactifi- 
cations can be found in Rcf. |3ql. 
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moduli space. Thus our results apply only at (or sufficiently near) those points in 
the moduli space. 



Heavy String Threshold Corrections 

We now discuss the heavy string threshold corrections in further detail. 

As we have discussed above, a natural unification of couplings occurs in string 
theory through which (at tree level in string perturbation theory) the gravitational 
(Newton) coupling constant Gn and all of the gauge couplings gi are related to one 
fundamental string coupling String- The precise relation is [ |31[ 

Q 

String = 87T ~^~ = 9ih foralH, (3.1) 

where a' is the Regge slope, and and k t are respectively the gauge coupling and 
Kac-Moody level of the gauge-group factor Gj. However, at the one-loop level, the 
above tree-level relations are modified to 

167T 2 167T 2 M s 2 triDg ~ 

— — = hi— + hi In ^ + Aj (3.2) 

9i [t 1 ) ^string f" 

where b{ are the one-loop beta-function coefficients, and Aj are the quantum correc- 
tions which reflect the contributions from the infinite tower of massive string states. 
These heavy string states are usually neglected in an analysis of the purely massless 
{i.e., observable) string spectrum, but they nevertheless contribute to the running of 
the gauge couplings. The threshold correction terms Aj therefore represent these con- 
tributions. Because these terms Aj enter these equations in the way that threshold 
corrections do in field theory, they are typically referred to as heavy string threshold 
corrections. 

There are several things to note about these threshold corrections. First, as we 
have already indicated, these threshold correction terms Aj are complicated functions 
of the moduli fields, and hence their values vary as one moves around in moduli 
space and considers different string vacua. However, our goal is to examine the 
role of these threshold corrections within a particular set of realistic free-fermionic 
models, and thus we are essentially evaluating these corrections at particular fixed 
points in moduli space. Consequently, for our purposes, these threshold corrections 
are merely numbers, and our conclusions concerning these numbers will apply only 
at (or sufficiently near) those points in moduli space. For our examination of gauge 
coupling unification within these string models, however, this will be sufficient. 

A second point concerns the renormalization scheme- dependence of these thresh- 
old corrections. Clearly the gauge couplings themselves are physical quantities, and 
thus their values completely scheme-independent. Consequently their unification is 
indeed a general property. By contrast, the values for the unification scale M string 
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and the threshold corrections Aj are scheme-dependent, and consequently their val- 
ues can only be specified within a particular scheme. For example, it is possible 
to choose a definition for M string which entirely absorbs the threshold dependence, 
leading to a definition of M stririg which depends on the VEV's of the string moduli. 
However, such a definition proves inconvenient, and for super symmetric theories it is 
preferable to work within the supersymmetric DR scheme. Thus, by selecting such a 
scheme, we have intrinsically selected not only certain definition for M string , but also 
a certain corresponding definition for the Aj which appear in the renormalization 



group equations ( |3. 2|) . We shall discuss these choices for the DR scheme below. 

Our final point concerns the dependence of these threshold corrections Aj on 
the gauge group (i.e., their dependence on the subscript i). In general, there are 
typically two terms that arise in a complete calculation of Aj, one of which depends 
on the gauge group in a complicated fashion (and which will be denoted Aj without 
the tilde), and the other of which depends on the gauge group only through the 
Kac-Moody level ki at which it is realized: 

Aj = A, + t t Y . (3.3) 



Here Y denotes the contribution which is independent of the gauge group [37| . How- 
ever, in our subsequent calculations, we shall only be interested in the relative running 
of the coupling constants g^, and this means that we shall only need to evaluate the 
differences between the threshold corrections Aj for different gauge group factors. We 
shall see this explicitly in Sect. 4. Moreover, since all of the realistic string models 
we will be examining have non-abelian gauge group factors realized at levels fcj = 1, 
the value of the group-independent quantity Y will be irrelevant for our purposes. 
We shall therefore need to evaluate only the simpler quantities Aj. 

Towards this end, the most important objects that we need to calculate are the 
one-loop string partition function Z(r) and the so-called "modified" one-loop string 
partition functions B g {t) corresponding to each gauge group factor G. Let us first 
focus on the partition function Z(t). In general, the partition function Z(t) of a 
given theory takes the form 

Z(r) = E("l) F Tr(«). (3.4) 

a 

Here the sum over a represents the sum over all sectors in the theory, the overall 
factor of (— 1) F ensures that spacetime bosonic and fermionic states contribute with 
opposite signs, and Tr(a) indicates a trace over the Fock space of mode excitations 
of the worldsheet fields: 

Tr(a) = Trg^g 17 " , q = e 2nir . (3.5) 

Here H a and H a are respectively the right- and left-moving Hamiltonians for the 
worldsheet degrees of freedom in the a-sector, and thus this trace simply counts 
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the number of string states at each worldsheet energy (H a ,H a ), as expected for a 
partition function. In string models this trace is generically realized as the result of 
a GSO projection between subsectors of the theory: 

Tr(«) = - c(a,f3)Tr(a,(3) . (3.6) 

9 p 

Here the /3-sum implements the GSO projection, c(a, (3) are the chosen GSO phases, g 
is a normalization factor, and Tr(a, (3) indicates a restricted trace over the appropriate 
(a, (3) subsector. 

The "modified" partition functions Bq which are needed in the calculation of 
the gauge coupling threshold corrections are then defined in a manner similar to the 
partition function, since they too must weigh the contributions from infinite towers of 
states. Indeed, they take the same form as Z in ( |3.4| ) except that they are multiplied 
by two factors of r 2 = Imr, and their corresponding traces are modified through 
insertions of the square of the spacetime helicity operator Q H and the square of the 
gauge group generator Qg- 

Tr(a) - Ti G (a) = r 2 2 Tr Q\ Q% q H «q^ (3.7) 

[with identical corresponding insertions into Tr (a, 13)]. Thus, these modified partition 
functions also count the numbers of states at each string energy level, but as expected 
each such degeneracy is multiplied by the gauge charge carried by the corresponding 
state. 

Finally, just as the one-loop vacuum energy is defined as the integral of the ordi- 
nary partition function over the modular-group fundamental domain T, the threshold 



correction contribution from the infinite tower of massive string states is given |19| 
by an analogous integral: 



A 



G 



? T 2 



B G (r) - r 2 b G 



(3.8) 



Note that since the contributions of the massless string states have already been 
included the calculation of the beta-function coefficient be appearing in (|3.2|), they 
must be subtracted from the modified partition function Bq so that A^ records 
only the extra contributions from the infinite towers of massive string states. This 
subtraction which excludes the massless states also renders the integral ( |3.8|) finite. 
Note that while the measure of integration in ( |3.8| ) is modular- invariant, the integrand 
is not. We shall discuss the modular properties of this expression in Sect. 6. 

As we have indicated at the beginning of this section, there are various assump- 
tions which enter the derivation of the one-loop result ( |3.8| ), and consequently there 
are various conditions under which (RTST) may be used. We shall merely list them 



here for completeness; a full discussion can be found in Ref. [19]. First, as mentioned 



above, this result is scheme-dependent, and has in fact been derived in the so-called 
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DR renormalization scheme; this is the scheme that is typically used for supersym- 
metric field theories and renormalizat ion-group equations. Therefore, as we shall see 
in Sect. 5, it will be necessary to include explicit scheme-conversion corrections when 
comparing with low-energy data obtained through other schemes (such as the usual 
MS scheme). Second, as we have discussed above, this choice of scheme in turn fixes 



the choice of string unification scale [19], and amounts to the definition 



e (l- 7 )/2 3-3/4 

^string = : ^string ^fpianck (3-9) 

47T 



where 7 ~ 0.577 is the Euler constant, and where String = y 327r/(a / Mp lanck ) is the 
string coupling. Numerically, this yields 



M string « £ string x 5 x 10 17 GeV . (3.10) 



Third, we reiterate that the result (|3.8|) is only a partial result which neglects the 
additional contributions such as the term Y in ( j3.3|) that are the same for all gauge 
group factors in a given model. Consequently, by using (|3.8j ) we cannot consider ab- 
solute values of the threshold corrections for a given gauge group factor, but rather 
only the relative differences of the threshold corrections Kq 1 — Ag 2 = — Aq 2 
between two gauge group factors G\ and G2 realized at equal Kac-Moody levels fcj.Q 
However, beyond these restrictions, the one-loop expression in (|3.8|) is completely 
general. In particular, it makes no additional assumptions about the structure of a 
given heterotic string model, or the values of various string moduli. Thus this expres- 
sion ( |3.8| ) will be our starting point when evaluating the threshold corrections in the 
realistic free-fermionic models. In particular, since these models have (2, 0) [rather 
than (2, 2)] worldsheet supersymmetry, we cannot use various moduli- dependent ex- 



pressions |33|, |34| which are derived from this. 



Threshold Corrections in Free-Fermionic Models 

We now concentrate on how this general procedure for calculating threshold cor- 
rections applies to the case of string models built through the free-fermionic con- 
struction || |10J. As briefly discussed in Sect. 2, in the free-fermionic construction of 
four-dimensional heterotic string models, the light-cone gauge worldsheet field con- 
tent consists of two transverse spacetime coordinate bosons, their two right-moving 
fermionic superpartners, and an additional set of 62 purely internal fermions of which 
18 are right-moving and 44 are left-moving. Collectively these 64 Majorana-Weyl 
fermions may be denoted tpi (£ = 1, 64), with an ordering such that £=1,2 corre- 
spond to the right-moving superpartner fermions carrying spacetime Lorentz indices, 

' We point out, however, that there has recently appeared an alternative procedure for calculating 
threshold corrections [Q which includes these gauge-independent constant terms. We will discuss 
this briefly at the end of Sect. 6. In any case, these constant terms will not be necessary for our 
analysis, so it will be sufficient for our purposes to use the definition (3 
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3 < t < 20 correspond to the purely- internal right-moving fermions, and t > 21 
correspond to the purely-internal left-moving fermions. A string model is then real- 
ized by specifying, for each sector of the theory, the boundary conditions of these 64 
fermions as they traverse the two non-contractible loops of the torus, as well as a set 
of phases which specify the generalized GSO projections which are to be applied in 
that sector. These parameters are not all independent, however, and must be chosen 
in such a way that certain self-consistency conditions (guaranteeing a proper world- 
sheet supercurrent, proper spacetime spin-statistics assignments, physically sensible 
projections, and modular invariance) are satisfied. The rules governing the construc- 
tion of such models are given in Refs. |5], [H], [TI|]. In particular, it is shown that 
the parameters describing the 64 fermionic boundary conditions in all sectors can 
ultimately be described via a simple "basis set" of (N + 1) different 64-component 
vectors V, (0 < i < N), and that the GSO phases in all sectors can be similarly 
described through a single matrix kij (0 < i,j < N). The number of basis vectors N 
which are necessary is of course mo del- dependent, increasing with the complexity of 
the model desired, and the complete constraint equations which relate self-consistent 
choices for the V,- vectors and ki 



Hj matrix are given in Ref. [11 
In the free- fermionic construction 



each of the a-sectors discussed above corre- 
sponds to a different set of fermionic worldsheet boundary conditions around the 
spacelike cycle of the torus, whereas the GSO projection is realized through the (3- 
summation over different fermionic boundary conditions around the timelike cycle of 
the torus. The corresponding GSO phases and traces are then unambiguously defined 
once a particular model {i.e., a set of Vj and ky parameters) is specified. For com- 
pleteness we now write these phases and traces explicitly in terms of the Vj and k 



i.r 



Note that henceforth we will be following the notation and conventions of Ref. [[LI 



where the definitions of any unexplained symbols below can be found. First, in the 
free-fermionic models, the traces Tr (a, (3) in ( |3.6|) are given by 



Tr (a, (3) = Tr [exp {-2mf3V ■ N' qV } q HaV q H ^ 

1 - u i 



20 



-1 



72 " =12^4 \ II \l 6 




(3.11) 



where N' aV indicates the worldsheet fermionic number operator in the aV sector with 
the Ramond zero modes excluded, the asterisk indicates complex conjugation, and 
where 



(if 



«V 



V (r) = g 1/24 n (1 - <f ) 

e 



n=0 
oo 



r) = V e 27Tibn q( n+a 1 2 / 2 



(3.12) 
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Note that the T2 1 factor in ( |3.11| ) represents the contribution to the trace from the 



zero-modes of the two transverse coordinate bosons X\ and that the sign of the 
square roots in ( ft.lip is not of practical importance for our calculation because there 



will always be an even number of v G factors for each type of O-function. This latter 
property is guaranteed by the so-called "cubic constraint" of Ref. [|TT[ . 

With the definition given in ( |3.11|) , the corresponding GSO phases and normal- 
ization factor in ( |3.6|) are then given by 



c(a 



(3) = r^exp^vrz^A^^ + ^-V.-aV)! 



N 



g = U m i ( 3 - 13 ) 



i=0 



where = ±1 represents the additional phase contribution which can arise due to 
permutations of the Ramond fermionic zero-mode gamma-matrices. The definition 
of is given explicitly! in Eq. (3.6) of Ref. [|TTJ. Note that = ±1, since for 



a given sector (aV, /TV) it can be shown that always contains an even number 
of factors of for each necessarily real fermion i. The quantities rrii in ( |3.13| ) are 



defined in Ref. 11 



The calculation of the "modified" partition functions Bq is similar to that de- 
scribed above for Z\ indeed, the GSO phases remain unchanged, and the only change 
is that the traces are modified as in Q3.7| ) due to the insertions of the charge opera- 
tors and two extra factors of T2. Now, these two extra factors of T2 are simple overall 
factors which merely increase by two the power of T2 in (|3.11| ). Let us therefore focus 
first on the gauge group generator Qq. In general, this object may be defined as 

Ql = ^ £ {Q {r) Y (3-14) 

where R is the rank of the group G and where are the normalized charge ele- 
ments of the Cartan subalgebra. In free-fermionic models, these charges Q gener- 
ally appear as linear combinations of the individual charge operators corresponding 
to individual worldsheet fermions; such linear combinations describe how the gauge 
group G is ultimately embedded or realized through free fermions. Note that in the 



* We point out, however, that when calculating T^-, one must include 7-matrix factors from 
only those Majorana-Weyl worldsheet fermions which are necessarily real, i.e., those fermions which 
cannot be globally paired in all sectors with any other Majorana-Weyl fermion in all sectors to form 
a complex Weyl fermion. This crucial restriction is not stated in the definition given in Ref. (T^j, 
but must be incorporated in the cases of models which simultaneously contain not only necessarily 
real fermions but also complex fermions with twisted "multi-periodic" boundary conditions. Most 
free-fermionic models of phcnomenological interest fall into this class. This restriction was also 
independently pointed out in Ref. |16f. 
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free-fermionic construction, the gauge group can be realized only through those world- 
sheet real Majorana-Weyl fermions which can be consistently paired in all sectors to 
form complex Weyl fermions, and it is indeed only for such pairs of real fermions 
- or equivalently for such complex Weyl fermions — that these fermionic charge 
operators can be defined. Let us write the indices of such complex Weyl fermions 
with capital letters, so that L = (£1,^2) represents a pair of real-fermion indices with 
corresponding charge operator Ql- We find, then, that one can generically expand 
Qq as a polynomial in terms of the underlying fermionic charge operators Ql'- 



Ql 



C LM QlQm 



(3.15) 



(G) 

where c L1 ^ are the model-specific coefficients which describe the gauge group embed- 
dings. 

Having related the gauge-group generator Q 2 G to the charge operators Ql oi the 
individual (complex) fermions, it is now straightforward to calculate the effect of 
the insertion of Q G into the trace. Let us first consider the effect of a single term 
Ql on the contribution to the trace from the corresponding L th complex fermion, 
where L = (A.) ^2)- Without the insertion, this isolated contribution from the two 
real fermions (£±,£2) would have been 



7/ 



-1 



'e 



+bt! 



e 



-ae 2 
+h 2 



rf x Q 



-a^ 
+&/, 



if 1 Yl exp {2mb h n} q {n - a ^ /2 . (3.16) 

n=— oo 



The first equality follows since the l\ and £2 fermions cannot be joined to form a 
complex fermion L unless they have the same boundary conditions in all sectors, so 



that 



a,£ 2 and b tl = bp 2 . The second line is just the expansion of the single 



function. Let us recall, however, that we may equivalently bosonize this complex 
fermion, in which case each term in this sum represents the contribution to the single- 
boson trace from a state in the corresponding bosonic one-dimensional momentum 
lattice P with lattice coordinate P = n — (so that the power of q in the expansion 
is the usual worldsheet energy contribution H = P 2 /2). But this bosonic momentum 
lattice P is nothing but the corresponding fermionic charge lattice Q ||. Thus we 
immediately recognize that the charge of each state contributing to the n th term in 
(|3.16|) is simply n — a^, so that the trace with a single Ql insertion is simply 



-1Q' 



-ai! 
-bp, 



1 Yl (n-a h )exp{2nib h n} q^^l 2 



(3.17) 



Likewise, an insertion of Q\ yields 



— 1 r\» 



—at! 
+b f , 



1 £ (n- a h ) 2 exp {2mb h n} q^- a ^ 2 l 2 



(3.18) 
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Thus for each term QlQm in Q3.15 ), the modified trace is the same as ( |3.11[ ) except 
that certain functions are replaced^ by either B' or 0" as defined above: for L = M 
we replace the corresponding single factor of with 0", and for L^Mwe replace 
each of the two corresponding factors of with 0'. 

— 2 

The helicity operator insertion Q H can be handled in precisely the same way, since 

the spacetime helicity operator is essentially the same as the gauge charge operator, 

only applied to the two real right-moving fermions L = (£1,^2) — (1)2) which carry 

spacetime Lorentz indices, and with an overall minus sign (to take into account that 

we are dealing with the corresponding complex conjugate ©-functions). Thus, the 
— 2 

insertion of Q H into the trace simply amounts to an analogous replacement of the 
corresponding functions with —0 , as discussed above. 

The final issue we discuss is a practical matter concerning the integration which is 
then necessary to obtain a value of Aq as in ( |3.8| ). It is of course possible to calculate 
be independently from explicit knowledge of the massless spectrum of the particular 
string model in question. However, there is an important self-consistency check which 
can be performed, since if we expand the corresponding modified partition function 
Bq as a power series of the form 

T 2 ECi?Y, (3-19) 



B, 



G 



(G) 

we see that the coefficient b 00 , which tallies the contributions from the massless 



states with total (left, right) worldsheet energies 
be equal to be- 



[n,m) = (0,0) respectively, should 

(3.20) 

Thus, the subtraction of T2&G within the integrand of ( |3.8| ) can be most easily achieved 



l(G) 



>G 



by expanding Bq in the form ( 3.19Q , and then simply setting the coefficient 1° 



(G) 

zero. In any case, as remarked earlier, a non-zero value of b 00 will result in a divergent 

(G) 

integral for A G , so setting the %, coefficient to zero is an efficient analytic way of 
removing what would otherwise be a logarithmic divergence. The resulting integral 
can then be performed by computing the sum 



a c = y b^ i^ 1 

<-r / j "ran m. 



) 

mn 



(3.21) 



§ In performing this replacement, there is in fact a subtle sign ambiguity: although 



—a 
b 



e 



as an algebraic identity, we find 0' 



-& 



a 

-b 



The question then arises as to 



whether 0' 



9' 



is to be used in the modified trace, amounting to an ambiguity in 



the overall sign of the charge Ql- However, after the GSO projections are performed, the resulting 
charge lattice Q must be always be invariant under the inversion Q — > Q (this is tantamount 
to CPT invariance). Thus, one must ensure only that consistent relative signs are used between 
GSO-related sectors. This amounts to choosing the same sign for ai in all sectors with identical 
a-boundary conditions but different /^-boundary conditions. 
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where the general integrals 



til = I ^Tq 11 . (3.22) 

can each be numerically evaluated. This procedure has the dual advantages of in- 
corporating the cancellations inherent in the GSO projections at an early stage (by 
calculating the coefficients prior to integration), and of analytically removing the 
divergence from the massless physical states % . Furthermore, the integrals (|3.22| ) 
then need to be evaluated only once, for any change in the particular string model 
under discussion merely amounts to a change in the corresponding coefficients b^. 
Note that inherent in our procedure is an exchange in the order of the summation 
over (m, n) and the integration over r. As discussed in Ref. [Q, this procedure is 
valid provided that the sum in ( |3.19|) is convergent as r — > ioo. However, this is 
necessarily the case for all string theories lacking charged physical tachyonic states. 
Indeed, any divergences arising from different sectors in the theory will have already 
cancelled in the sum over all sectors, so that the total coefficients b mn which we use 
necessarily have b nn = for all n < 0. Hence, once the contributions from all sectors 
have been added together (as we have done here), the interchange is valid. 

Finally, we point out that in removing the logarithmic divergence from the mass- 
less states, we have defined &oo through the expansion ( [3.19| ). In particular, we have 
not followed Ref. |T9[ in defining 



*ioc 



lim \ t 2 l B(t)\ , (3.23) 



as there can be cases for which such an identification is incorrect. For example, if 
B(t) has non-zero coefficients b mn with m+n < and m ^ n {i.e., contributions from 
charged "unphysical tachyons"), then the r — ► zoo limit of B(t) diverges, and ( |3.23p is 



in error. We shall see explicit examples of this sort of occurrence in Sect. 4, where we 
shall be considering certain models which are free of physical charged tachyons (i.e., 
which have b nn = for all n < 0), but which nevertheless have non-zero coefficients 
b mn for negative unequal values of m and n. Although we shall find that m + n > 
for such coefficients, the appearance of such contributions with negative and unequal 
values of m and n shows that these sorts of unphysical tachyons generically appear, 
and that unphysical tachyons with m + n < cannot be generically ruled out. In 
any case, however, such unphysical tachyons do not lead to the sorts of divergences 
in the threshold corrections A that would arise from physical tachyons with m — n, 
since the contributions from unphysical tachyons are projected out of the integral in 
the r — > oo region. We shall discuss the implications of such unphysical tachyons in 
Sect. 6. 
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4 Threshold Corrections in Particular Models 



In the previous section we discussed how threshold corrections A^ can be calcu- 
lated in general free-fermionic models. In this section we now apply this procedure 
to the cases of the realistic free-fermionic models described in Sect. 2. In particular, 
we shall explicitly calculate the threshold corrections A^ for each gauge group factor 
in each realistic free-fermionic model which has appeared in the literature to date. 
This includes the "revamped" flipped SU(b) model, several SU(3) x 577(2) x U(l) 
models, and an £0(6) x 5*0(4) model. Furthermore, to test the dependence of our 
results on the existence of spacetime supersymmetry in these models, we shall also 
consider the case of a certain ncm-supersymmetric SU(3) x 577(2) x U(l) model. 

As we discussed in the previous section, our procedure for each model is essentially 
the same, and can be broken down into several distinct steps. 

First, for each model, we must determine the embedding of each of its gauge group 
factors within the free worldsheet fermions — i.e., we must essentially determine the 
coefficients which appear in (|3.15|) . We will refer to this procedure as obtaining 
the appropriate charge polynomial for each factor. In general, this requires not only 
analyzing how each gauge group factor arises within the free-fermionic construction 
for each model, but also calculating the appropriate overall normalizations. These 
issues will be discussed for each model below. 

Second, for each model, we must then actually perform the threshold calculations, 
and calculate the sum over all of the contributing sectors that appear in the model. 
For each of the realistic models that we will be analyzing, this typically amounts to 
summing the separate contributions from several thousand individual sectors for each 
gauge group factor. 

Third, for each model, we then perform a set of important self-consistency checks 
by verifying that the relation ( |3.2U| ) holds for each gauge group factor in the model. 
In order to verify this relation, we therefore have to do two separate calculations: 
we must calculate the 6qo coefficient that arises from an explicit expansion of the 
modified partition function B^ g \t) for each gauge group factor G, and we must also 
independently calculate the beta-function coefficient be using our prior knowledge of 
the massless spectrum of the model in question. Ensuring that these two quantities 
agree for all gauge group factors thus guarantees that our calculations of both the 
massless spectrum of the model and the threshold corrections from its massive states 
are not in error (or inconsistent with each other). This also provides a non-trivial 
check of our charge polynomials, which in turn verifies our determinations of both 
the gauge group embeddings, and their non-trivial normalizations. 

Finally, for each model, we then put all of these results together to calculate the 
resulting threshold corrections Aj. The effects of these threshold corrections on gauge 
coupling unification will be analyzed in Sect. 5. 
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4.1 Flipped SU(5) model 



We first analyze the threshold correction in the "revamped" flipped SU(5) model 
pOj] . The defining parameters (worldsheet fermion boundary conditions and GSO 
projection phases) which yield this model are listed in the Appendix. The analysis 
of heavy threshold corrections for this model was already performed in Ref . |39| , but 
we shall use this example to set the procedure for the analysis in subsequent models. 
Furthermore, as we shall see, our results differ substantially from the numerical results 



of Ref. |39| . However, we believe that this difference is ultimately a reflection of 
our improved method of calculating the modified partition functions B(t), and of 
numerically integrating these functions B(t) over the modular-group fundamental 
domain. In particular, as discussed in the previous section, our procedure analytically 
removes the logarithmic divergence from the massless states, and in so doing explicitly 
verifies that be and 6qo agree (as required for self-consistency). Furthermore, we will 
see that the size of our results is more in line with those from the other realistic free- 
fermionic string models we will be examining, as well as from previous string threshold 



calculations in various orbifold |L9], [40[ and Type-II |^T| models. We therefore believe 



that our results are more reliable. 

Determining the Charge Polynomials 

This model has the observable gauge symmetry SU(5) x U(l), which is realized 
through the following five complex (left-moving) worldsheet fermions: 

L — l ^ i t 
L = 2 l t 

L = 3 « l t 
L = A 

L = 5 li = (37,53) < — fermion ^ . (4.1) 

Explicitly, this means that we are labelling as L = 1, 5 the five complex fermions 
which can be formed, as indicated, from the real fermions 33, ...,37 and 49, 53. This 
latter numbering reflects the ordering of the real fermions as listed in the Appendix. 
These five complex fermions L = 1, 5 are respectively the complex fermions called 
ip in Sect. 2. 



(33,49) «- 


— fermion ip l 


(34,50) «- 


— fermion ip 2 


(35,51) «- 


3 

— fermion i/j 


(36,52) <- 


— fermion ip 4 


(37,53) «- 


— 5 

— fermion ijj 



As discussed in Ref. |39|, the four traceless £77(5) generators for this model can 
then be written in terms of the charges Q l corresponding to these complex fermions, 



2.3 



Q^SU(5) ~ ~TTK (^1 + Ql + Qz — 3Q4) 



12 

Q§ (5) = _L(Q 1 + Q 2 + Q 3 + Q 4 -4Q 5 ) , (4.2) 



while the remaining single orthogonal U(l) generator (essentially the trace of the 
original larger £7(5) symmetry) is therefore 

Q§ x) = -±=(Q l + Q 2 + Q 3 + Q A + Q 5 ) . (4.3) 

Note that each of these generators is normalized so as to produce conformal dimension 
one for the massless states. 

Given these individual generators for the SU(5) and U(l) group factors, we 
then compute the corresponding squared polynomials Q 2 G for each factor according to 
(ft. 14 ). Note that since the complex fermions L = 1, 5 all have identical boundary 



conditions in all sectors of this model, their charges Ql will be identical in all sectors. 
It is therefore possible to replace any occurrence of products of charges of the form 
QlQm with Q 2 if L = M, and with Q1Q2 ii L M.f\ The resulting squared 
polynomials then take the simple form 

QsU(5) = Ql ~ Q1Q2 

Qu W = Ql + 4QiQ2 ■ (4.4) 

Thus, the only two "basis" charge insertions that we need to consider are Q± 2 and 
Q1Q2. 

Note that since we will ultimately be interested in only the difference of the 
threshold corrections for SU(5) and U(l), we could at this stage proceed to consider 
only the difference of their corresponding charge polynomials in (|4.4j ). Indeed, this 
would cancel the Qi 2 term completely, and leave us with a single "basis" insertion 
Q\Q2- However, this simplification would later rob us of an important self-consistency 
check on the corresponding beta-function coefficients, for we would be able to compare 
only the difference of these coefficients, rather than each one separately. For this 
reason we shall retain both charge polynomials in our basis. 

Threshold Contributions: Results of Calculation 

Following the procedure outlined in the previous section, we now calculate the 
separate contributions to the threshold corrections that arise from each of these "ba- 
sis" insertions. Specifically, for each basis insertion QlQm (i-e., either Qi 2 or Q1Q2), 
we calculate the corresponding trace B{r) over mass levels and all sectors, its massless 

* Note, however, that we must continue to distinguish between the squared charge for a single 
fermion, and the product of single charges for two separate fermi ons. These two charge combinations 
have different effects within the trace, as discussed below ( [3.1g ) . 
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expansion coefficient 600 > and its corresponding integral A with its massless coefficient 
600 set to zero (to analytically remove the divergence from the massless states). Our 
results follows: 

' b 00 A 



~Q7~ 

Q1Q2 



3.5 12.331 (4.5) 
4.5 1.53625 



In evaluating A, we have followed the procedure indicated in (|3.21| ) and (|3.22|) , 



analytically g-expanding each modified partition function B(r) to the fourth excited 
energy level, and numerically evaluating the convergent integrals I$ n to the precision 
indicated in the results quoted above. This procedure thus minimizes any numerical 
error, so that essentially none remains to this level of accuracy. 

Self- Consistency Checks 

Next we perform our self-consistency checks by comparing these values of 600 for 
each charge insertion with the one-loop beta-function coefficients expected from the 
(known) spectrum of this model. In general, the one-loop beta-function coefficients 
are given by 

bsu(5) = -15 + 2N g + N h 
b u{1) = 2N g + N h + ^N 4 , (4.6) 

where N g is the number of generations and antigenerations [with a generation con- 
sisting of a full 16 representation of SO(10) decomposed under SU(5) x U(l)], where 
Nh is the number of 5 + 5 Higgs pairs, where iV 4 is the number of hidden SU(4) pairs 
4 + 4 that carry U(l) charge ±5/4. In the revamped flipped SU(5) model, there are 
four 16 representations and one 16 representation of 5*0(10), Nh = 4 and N& = 6. 
Thus, the expected one-loop beta-function coefficients are 

bsu® = -1 , b u{x) = 21.5. (4.7) 



We now compare these numbers with the values of 600 m (EH)- Using the charge 
polynomials given in ( |4.4| ), we indeed find 

bsu(5) = 3.5 - 4.5 = - 1 
b m = 3.5 + 4(4.5) = 21.5 , (4.8) 

in perfect agreement with (|4.7|). 



Final Threshold Corrections 

Given these results, we now compute the final relative threshold correction Aym — 
As[/(5). Since - Q 2 SU ^ = 5QiQ 2 , the difference Au(i) - ^su(5) will simply be 

five times the value of A in ( |4.5| ) for Q1Q2, or 

A £ /(i) - Asu(s) = 7.68125. (4.9) 
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Note that this result is smaller by approximately a factor of three from the value that 



was found in Ref. p9 |. 



Given this relative threshold correction, we can then compute its effect on the 
string unification scale M string . From the renormalization-group equations ( |3.2|) for 
the SU(5) and U(l) couplings, and by taking the difference between the 577(5) and 
£7(1) one-loop RGE's, we find that the corrected unification scale is in general given 

by ' 

M£r ed) = ^~ ted) exp|^^} • (4.10) 
Since bum — bsu(5) — 22.5, we find using ( |4.9|) that 

string 



^s ( tring CtCd) ~ 5.93 x 10 17 GeV (4.11) 



for this model. 

We remind the reader that in the realistic free-fermionic models, one combination 
of the £7(1) factors is anomalous. This anomalous £7(1) combination arises due to the 
fact that these free-fermionic models are (2,0) compactifications rather than (2,2). 
The anomalous £7(1) gives rise to a Fayet-Iliopoulos D-term that breaks supersymme- 
try and destabilizes the vacuum. The models must therefore choose non-zero VEV's 
for some of the scalar fields (twisted moduli) so as to cancel the anomalous £7(1) 



.D-term | 42|| . This corresponds to a shift in the string vacuum, and one can then ask 



whether this shift can modify the estimate of 8 A, and consequently our estimate of 
Muring • However, the contributions to the threshold corrections come from string 
states weighing > g s tring^pianck/v / 8/T, whereas the (extra) masses acquired from shift- 
ing the string vacuum are of higher order in String- Hence, these extra masses only 
affect the value of Mg t c ° n r g e€ted ^ at higher order, and can be ignored in our analysis [[3fJ . 
The same also holds true for the other models we will be examining. 

4.2 First SU(3) x SU(2) x U(l) model 

This model [[H]] is listed in the Appendix. The observable gauge group of this 
model is SU(3) x SU(2) x £7(1) 2 , with one of the £7(1) factors associated with the 
SU (3) co i OI group factor (arising as the trace of the original larger Z7(3) color sym- 
metry), and the other associated with the S , [/(2) le f t group factor (and arising as the 
trace of the original larger £7(2) weak isospin symmetry). We shall denote these two 
£7(1) factors as U(l)c and U(1)l respectively. The electroweak hypercharge £7(1) is 
then a linear combination of these two £7(1) factors. 

Determining the Charge Polynomials 

All of these group factors are realized through the same five complex fermions as 
in Q4.1D for the flipped SU(5) model. Unlike the case of the flipped SU(5) model, 
however, these five fermions no longer share the same boundary conditions (indeed, 
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it is for this reason that the gauge group is altered). Rather, only the L = {1, 2, 3} 
fermions share the same boundary conditions; likewise, the remaining L = {4, 5} 
fermions share the same set of different boundary conditions. Therefore, we expect 
that our needed charge polynomials must be built from the larger set of five elemen- 
tary "basis" charge bilinears {Qi 2 , Q1Q2, Q1Q4, Qi 2 , QiQs}- 

Explicitly, the generators of the SU(3) x SU(2) x U(l) 2 group are as follows. The 
traceless £77(3) generators are given by 



1 



Qsu { 3) = ^ (Qi + Q2 ~ 2Q 3 ) , (4.12) 

and the orthogonal U(l) c generator corresponding to this color group factor is the 
trace of the original U(3) c symmetry: 

Qu\i) c = -^(Qi + Q2 + Qs) • (4.13) 

As before, we are properly normalizing all of the linear combinations of charges so 
that massless states will have conformal dimension one. Since the L — 1,2,3 fermions 
all have identical boundary conditions in every sector, the squared charge polynomials 
corresponding to these factors are then simply given by 

Q 2 su(3) = Qi ~ Q1Q2 

Q 2 u(Dc = Qi + 2 QiQi ■ ( 4 - 14 ) 

Likewise, the traceless SU (2) generator is 

Qsu { 2) = ^(Q 4 -g 5 ) , (4.15) 

and the corresponding orthogonal U(1)l generator is the trace of the original U(2)l 
symmetry: 

Qu(i) L = ^(Q4 + Q 5 ) • (4.16) 

Because the L = 4, 5 fermions also have identical boundary conditions in all sectors, 
the polynomial insertions in the modified partition function are simply given by 

Qsif(2) = Qi ~ Q4Q5 

In this model, the electroweak hypercharge is a combination of U(l)c and U(1)l, 



QliDr = Ql + • (4.17) 



U(1)y = \u{l)c + \U(\) L , (4.18) 
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where the coefficients in this linear combination are appropriate for the un-normalized 
traces U(l)c and U(1)l- Thus, the properly normalized charge polynomial corre- 
sponding to U(l)y is simply 

6 ( Q1 + Q2 + Q3 Q ± ±Q5\ 2 
5 V 3 2 J 

\(^Q\ + \QiQ2 + 2QiQ* + \qI+\q4Qs) , (4.19) 

where in the second line we have again used the fact that the boundary conditions 
within each of the fermion sets L = 1, 2, 3 and L = 4, 5 are identical. 

Threshold Contributions: Results of Calculation 

Given the charge polynomials Qs U3 \, Qsu(2)i an d Qum Y f° r this model, we see 
that, as anticipated, we need to calculate the contributions from all five of the "basis" 
charge insertions {Qi 2 , Q1Q2, Q1Q4, Q4 2 , QiQb}- Following the same procedure as 
outlined in the previous section, we calculate the 600 coefficients and integrals A 
corresponding to each insertion, with the results: 





boo 


A 


Qi* 


3.5 


9.51642 


Q1Q2 


4.5 


1.53625 


QiQa 


0.5 


-0.291404 




9.5 


12.951 


Q4Q5 


4.5 


1.53625 



(4.20) 



As an aside, we note that the modified partition function B for the insertion of 
Q1Q2 is the same as that with the insertion of This explains their identical 

values of 600 and A. Moreover, this is also the same function B which appeared for 
the Q1Q2 insertion in the flipped SU(5) model. This is a reflection of the similarity 
of their underlying free-fermionic structures. 

Self- Consistency Checks 

Let us now verify the self-consistency of the above values of b 00 . Indeed, as we 
shall see, these checks become increasingly non-trivial as the models become more 
complex. 

For this model, the full massless spectrum was presented in Ref. ||21|| . Here we 
list only the states and their non-trivial charges under SU(3)c x SU(2) L x Z7 (l)c? x 
U(1)l- The sectors b 1; b 2 , and b 3 produce three 16 representations of SO(10), 
decomposed under SU(3)c x SU(2)l X U(l)c x U(1)l- The Neveu-Schwarz sector 
produces three pairs of electroweak doublets with (U(1)c,U(1)l) charges (0, ±1). 
The sector bi + b 2 + a + [3 produces an additional pair of electroweak doublets and 
a pair of color triplets with (U(l)c, U(1)l) charges (0,±1) and (±1,0) respectively. 



Qu(i) Y ~ 
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Denoting the sector 1 + k>i + b 2 + b3 as J, and introducing the notation (+/) to 
indicate the two sectors with and without the vector / added, we find that the four 
sectors b 2 + b 3 + a ± produce an additional pair of color triplets with U{1) 

charges (±|,±|), and eight SU(3)c x SU{2) L x U(l)y singlets with U(l) charges 
(±|, =f|)- The sectors bi + b 3 + a ± 7(+/) produce an additional pair of electroweak 
doublets with £7(1) charges (±~,±~), and fourteen SU(3)c X SU{2) L singlets with 
U(l) charges (±|, =f|). Finally, the sectors bi + b 23 + a ± 7, b 2 + b 3 + a ± 7, and 
bi + b 2 + b 3 + a + (3 ± 7(+J) produce a total of 32 S77(3) c x SU(2) L singlets with 
U(l) charges (±f,±§). 

Given this matter content, the one- loop beta-functions coefficients are determined 
as follows. In general, we have 

bsu(3) = -9 + 2N g + N 3 
bsu(2) = —6 + 2N g + N 2 

bu(i) Y = JT^QuiDy (4.21) 

where N g is the number of 16 representations of S'O(IO) in the massless spectrum, iV 3 
in the number of color triplets in vector-like representations, and N 2 is the number 
of electroweak doublet pairs. The trace in the U(l) beta-function is taken over 
the entire massless spectrum, and the U(1)y normalization is fixed by the standard 
SO(10) embedding. Since U{l)c and U{1)l have the standard SO(10) embedding, 
U(1)y has the standard 5*0(10) normalization, ki = 5/3. Using N g = 3, iV 3 = 2, 
N 2 = 5, and the U(l) charges given above, we then find 

bsu(3) = -1 ; b S u(2) = 5 , 

b u(1)c = 12.5 , b u{1)L = 14 , b u[1)Y = 14.6 . (4.22) 

It is now straightforward to compare these five results with the values of 600 
listed in ( f4.20| ). Indeed, given the expressions for the charge polynomials in (|4.14j ), 
(|4.17|) , and ( |4.19|) , we find that in each case, the appropriate linear combinations 



of these values of 600 agree with the values obtained from the massless spectrum. 
This check is extremely non-trivial, essentially verifying not only the known massless 
spectrum and our evaluation of the modified partition functions B(t) as sums over 
all of the (thousands of) sectors, but also verifying the U(l) normalizations which we 
determined through other means. This is therefore an important consistency check 
on our analysis. 

Final Threshold Corrections 



Given the charge polynomials in ( [4.14| ), ( f4.17| ), and ( |4.19| ), we take the appropriate 



linear combinations of the five values of A listed in ( [4.20|) in order to obtain the 
relative string threshold corrections for the SU (3)c, SU (2) l, and U(1) Y group factors 
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(where Y refers to the normalized weak hypercharge). In this way we obtain the 
relative thresholds 

Ac (% - A 8U{3) = 5.0483 , A u{1) . - A su(2) = 1.6137 . (4.23) 

In Sect. 5 we will examine the effects of these threshold corrections on the experi- 
mental parameters. 

4.3 Second SU(3) x SU(2) x U(l) model 

We now turn to a different 377(3) x SU{2) x U(l) model |T7|, one whose 
analysis is substantially more complex due to the fact that portions of its gauge 
group are realized as enhanced symmetries. In particular, this means that some of 
the corresponding gauge bosons originate not in the Neveu-Schwarz sector of the 
theory, but rather in various additional twisted sectors. Therefore, for the sake of 
clarity, we shall first discuss the origin and structure of the gauge symmetry in this 
model before proceeding with the analysis. The parameters defining this model are 
given in the Appendix. 

The Gauge Structure of the Model 

Due to the enhanced gauge symmetry of this model, it turns out that the gauge 
group is ultimately realized not only through the five complex fermions listed in (|4.1|) , 
but also through the following additional five fermions denoted L = 6, 10: 



7 = 6 ti = (25, 31) < — fermion y 3 + vf 

7 = 7 <^=> li = (21, 30) < — fermion y 4 + iZJ 5 

7 = 8 h = (24, 28) < — fermion u 2 + iuJ 4 

7 = 9 it = (41, 57) < — fermion 1 

L = 10 ii = (48, 64) < — fermion 8 . (4.24) 



As before, the indices ii refer to the ordering of the real fermions presented in the 
Appendix, and the labels y 1 , ZU\ and 0* refer to the discussion in Sect. 2. Three other 
complex fermions which we will also need to consider in our discussion (although not 
in our eventual calculations) are 

L = 11 £i = (38,54) < — fermion rf 

7 = 12 ti = (39, 55) < — fermion rf 

L = 13 ti = (40, 56) * — fermion rf . (4.25) 

We begin by outlining the various U (1) factors which will appear. First, as before, 
there are the U(l)c and U(1)l factors which are respectively the traces of the larger 
U(3)c and U{2) L symmetries. These correspond respectively to the L = {1, 2, 3} and 
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L = {4,5} sets of complex fermions. Next, there are £7(1) factors which correspond 
to each of these additional complex fermions L = 6, ...,13 individually. We shall 
follow the literature in naming these £7(1) 's as follows: 



(L = 6) 
(L = 7) 
(L = 8) 
(L = 9) 



17(1)4 , 
U(l) 5 , 



(L = 10) 

(7 = 11) 
(L = 12) 
(L = 13) 



U(l) 9 , 

U(l)i , 

U(lh , 

U(l) a • 



(4.26) 



In this model, the observable gauge group formed by the gauge bosons from the 
Neveu-Schwarz sector alone is 



SU(3) C x SU(2) L x 17(1)0 x U(1) L x £7(1) ^,3,4,5,6 



(4.27) 



However, in this model a new feature arises due to the appearance of two additional 



gauge bosons from the twisted sector 1 + a + 27 p3| , and two corresponding new 
generators. These new gauge bosons are singlets of the non-Abelian group, but carry 
£7(1) charges. However, referring to these new generators as T , it turns out that we 
can define the linear combination 



rpo 



U(l) c + U(l) 4 + £7(1)5 + £7(l)e + £7(1) 7 - £7(l)g 



(4.28) 



in such a way that the three generators {T S ,T^} together form the enhanced sym- 
metry group SU (2). Due to the custodial role played by this additional SU (2) in this 
model, we shall refer to this new factor as SU(2) cust . Thus, we find that the original 



observable symmetry group ( 4.27 ) of this model has been enhanced to 



SU(3) C x S£7(2) L x S£7(2) cust x £7(l) c , x £7(1) L x £7(1) 1)2)3 x £7(1) 



4',5',7" 



(4.29) 



where (again following the nomenclature in the literature) we have chosen to de- 
fine the following linear combinations of remaining £7(1) factors, each of which is 
orthogonal to T 3 : 



U(l) c < 

£7(1) 4 ' 
£7(1) 5 ' 

£7(1) 7 „ 



\u{l)c-\u{l) 7 + \u{l), 

£7(1) 4 - £7(l)s 
£7(1) 4 + £7(1) 5 -2£7(1) 6 
5 



c 



£7(1) 4 + £7(1) 5 + £7(l) f 



+ £7(1) 7 -£7(1) 9 . (4.30) 



The final issue is the definition of the electroweak hypercharge in this model. In 
fact, due to the extended symmetry, we now have the freedom to define the weak 
hypercharge in several ways. One option is to define the weak hypercharge just as in 
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( |4.18| ) for the previous SU(3) x SU(2) x U(l) model. However, in the present model, 
the U(l)c symmetry is now only part of the extended custodial symmetry SU (2) cust . 
Indeed, expressing U(l)c in terms of the new linear combinations defined above, we 
have 



T A + - U r > 
5 



(4.31) 



Thus U(1)y, by depending on the T 3 of the custodial SU(2), is no longer orthogonal 
to SU(2) cnst . We must therefore instead define the new linear combination with this 
term removed, 



U(1) Y , = U(l), 



1 



1 



f/(l) 4 + f/(l) 5 + [/(l) 6 + [/(l) 7 -f/(l) £ 



so that the weak hypercharge is expressed in terms of U(1)y' as 

1 



U(1) Y = U(1) Y , + -T 3 



Qe 



T'l + Y = Tl + Y' + - T 3 



1 

7 

2 cust 



(4.32) 



(4.33) 



The final observable gauge group then takes the form 



SU(3) C x SU(2) L x SU(2) cust x U(1) Y > x j seven other 17(1) factors j . (4.34) 

Of course, these remaining seven U(l) factors must be chosen as linear combinations 
of the previous U(l) factors so as to be orthogonal to the each of the other factors 
in ( PI . 



Determining the Charge Polynomials 

Given the gauge structure outlined above, we now determine the corresponding 
charge polynomials for the physically relevant SU (3)c, SU (2) L , SU(2) cust , and U(l) Y i 
factors. 

The charge polynomials for the color SU(3)c and the electroweak SU{2)l gauge 
groups are the same as in the previous SU(3) x SU(2) x U(l) model: 



Q 2 



Q 2 



SU(3) 
2 

5(7(2) 



Qi ~ Q1Q2 
Q 2 — Q4Q5 ■ 



(4.35) 



The charge polynomial for the custodial SU (2) cust gauge group is obtained from the 
expression (|4.28|) for its diagonal generator, and, with proper normalization, is given 
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by 
Q 



2 

SU (2) cus t 



- \3Qi + 6 (Q 1 Q 2 + Q X Q & + Q 1 Q 7 + Q X Q % + Q ± Q 9 - Q 1 Q 10 ) 



+ Ql + 2 (Q 6 Q 7 + Q 6 Q S + Q 6 Q 9 - Q 6 Q 



io 



+ q? + 2 (Q 7 Q 8 + Q 7 Q 9 - Q 7 Q 

+ — 2Q9Q10 + Qioj • 



2 (QsQb ~ QsQio) 
(4.36) 



Likewise, the normalized charge polynomial combination for U(1)y> is found to be 
Ql 



fu(i) Y , 



24 f 2S S 
T ^{ T ^(Qi + 2Q 1 Q 2 ) + -Q 1 Q, 



- —(QiQa + Q1Q7 + QiQs + Q1Q9 - Q1Q10) 
1 

+ 64 



Ql + 2(Q 6 Q 7 + Q 6 Q 8 + Q 6 Q 9 - Q 7 Q 10 ) + Q 2 7 
+ 2(Q 7 Q 8 + Q 7 Q 9 - Q 7 Q 10 ) + Qj 
+ 2(Q S Q 9 - Q S Q 10 ) +Ql- 2Q 9 Q 10 + Q 10 j . 



(4.37) 



In assembling these charge polynomials we have made use of identical boundary 
conditions wherever possible in order to simplify these expressions. 

Threshold Contributions: Results of Calculation 

Given the above polynomials, we see that we must now calculate the contributions 
from a basis set of 30 different charge insertions corresponding to the ten complex 
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worldsheet fermions L 



10. Our results are as follows: 
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(4.38) 



Note that unlike the simpler previous cases, some of these charge insertions re- 
sulted in modified partition functions B(t) containing contributions from charged 
unphysical tachyons. Specifically, in the g-expansion of some of these functions B(t), 
various coefficients b mn with m < or n < are non-zero. Those cases are indicated 
with an asterisk following the corresponding value of &oo- This occurrence is not 
unexpected, however, since such unphysical tachyons are generically present (and in 
fact unavoidable) in generic string models, and are required for the consistency of 
the theory ]4"5|] . Indeed, they do not lead to any divergence in the corresponding 
threshold integrals A, since the fact that they are unphysical (i.e., with m 7^ n, 
m — n G Z) implies that they have no contributions from the region 72 > 1 of the 
fundamental domain T from which infrared divergences might arise. We point out, 
however, that any unphysical tachyonic contributions with m + n < will render 
incorrect the expression ( |3.23D which was proposed in Ref. ||19|| . We will discuss some 
further implications of these unphysical tachyons in Sect. 6. 

Also note from the above results that the charges Qq and Q7 appear interchange- 
able insofar as their effects on the modified partition functions B(t) are concerned. 
Despite this fact, however, it can easily be verified that the corresponding L = 6 and 
L = 7 worldsheet fermions do not share the same boundary conditions. 

Finally, we observe that in the cases of the insertions Q^Qq, QaQt, and Q^Qs, 
the corresponding modified partition functions B(t) actually vanish identically. This 

occurs because each term in these expressions B(r) contains a factor of 0' ^' 







e' 





1/2 



or 6' 



1/2 




As can be seen from their definition in ( [3.17D , these particular 



singly primed B'-functions vanish. 
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Self- Consistency Checks 

Given the complexity of this model, it is crucial now more than ever to verify 
that our self-consistency checks are satisfied. The complete massless spectrum of the 
model is given in Ref. P3| . The model contains a total of eighteen color triplets, 
twelve from the sectors bi, b2, and b3 (which produce the light generations), and 
an additional three pairs from the sectors bi^ + b3 + (3 ± 7 and 1 + a + 27. Thus, 
from (|4.21|) , we see that the one-loop beta-function coefficient b S u{3) c for the color 
SU (3) gauge group factor vanishes. Comparing this against the above values of b 00 
(in particular, the difference between those for the Q\ and the Q1Q2 insertions), we 
see that agreement is obtained. 



Turning now to the SU(2)l electroweak symmetry, we recall [23[] that this model 
contains 28 electroweak doublets. Five doublets are obtained from each of the sectors 
bi, b 2 , and b 3 . An additional doublet beyond the Minimal Supersymmetric Standard 
Model arises due to the transformation of the lepton left-handed doublet under the 
custodial SU(2) cust symmetry. The Neveu-Schwarz sector and the sector bi + b 2 + 
a+/3 contribute three and two additional pairs, respectively. Finally, three additional 
electroweak doublets are obtained from each sector, 1 + bj + a + 27 for i = 1, 2, 3. 
Thus, from ( |4.21|) , we find b S u{2) L = 8. Consulting the above charge polynomial 
expression for SU (2) L and the corresponding entries in ( |4.38[ ), we see that once again 
the values agree. 

Finally, we compare the one-loop beta-function coefficients for the custodial 
SU (2) cust symmetry. It turns out that this model contains twelve massless doublets 
under the custodial S77(2) cust symmetry. Thus, the one-loop beta-function coeffi- 
cient of the SU(2) cnst group vanishes. From the polynomial expression ( f4.36| ) and 



the appropriate entries from fl4.38|) , it can be verified that the appropriate linear 
combination of 600- values vanishes as well. 

Final Threshold Corrections 

First we must discuss the Kac-Moody factors associated with the U(l) factors in 
this model. In the class of string models we have been examining, the Kac-Moody 
level of the non-Abelian group factor is always one. The situation is somewhat more 
complicated for the U(l) factors, however. In general, a given U(l) current U will 
be a combination of the simple worldsheet U(l) currents Uf = /*/ corresponding 
to individual worldsheet fermions /, and will take the form U = Y^f a fUf where the 
a,f are certain model-specific coefficients. The Uf are each individually normalized 
to one, so that (Uf,Uf) = 1. To produce the correct conformal dimension for the 
massless states, each of the U(l) linear combinations U must also be normalized to 
one. The proper normalization coefficient for the linear combination U is thus given 
by N = (J2f a /) _1//2 , so that the properly normalized U(l) current U is given by 
U = N-U. 

Now, in general the Kac-Moody level of the U(l)y generator can be deduced from 
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the OPE's between two of the U(l) currents, and will be 



k x = 2N~ 2 = 2 a} . (4.39) 

/ 

For a weak hypercharge that is a combination of several £7(1) 's with different nor- 
malizations, the result ( f4.39|) generalizes to 

h = ( 4 - 4 °) 

i 

where the ki are the individual normalizations for each of the f/(l)'s. 

Now, in the model analyzed in Sect. 4.2, the U(1)y generator is given as a com- 
bination of simple worldsheet currents that produces the correct weak hypercharges 
for the Standard Model particles. Thus, in that case, k\ is simply given by fl4.39| ). 
For the weak hypercharges ( |4.32j ) and (^4.33|) that appear in this model, however, we 
instead use ( |4.4U| ). Thus, for this weak hypercharge, we see from ( |4.33| ) and ( |4.40| ) 
that k\ = (l/4)/c 2c + ky = 1/4 + 17/12 = 5/3, which is the same as the standard 
5*0(10) normalization. 

From ( [4.14| ) and ( |4.17| ), we then obtain the following relative string heavy thresh- 
old corrections for SU(3)c, SU(2)l, SU(2) cust , and £/(l)y,: 

A SU (2 )l - A su{3)c = 2.2690 , A m . - A su{3)c = 5.3767 

A u{lh) - A su{2)cust = 2.8587 , A u{1) , - A su(2)cust = 2.4299 . (4.41) 

Note that for future convenience [see ( |5.26| )], we have defined 

Ac/( % = I (\ A su(2 )c + ^u(i hl ) , (4.42) 

whereupon 

17 

A u{1)jr - A SU (2) cust = — (A u{1)i , f - A SU (2) cust ) ■ (4.43) 
Alternatively, we can define the weak hypercharge to be the combination 

U(1) Y = \u{1)l + U{l) c > (4.44) 

where U(l)c' is given in (|4.30|) . In this case the Kac-Moody levels of U(1)l and 
U{\)c> are 4 and 5/3 respectively, so that k\ = 8/3. In this case we find that Asu(3) c 
and Asu(2) L are the same as above, and that the relative threshold for the properly 
normalized U(l)y is 

A u{1) . - A su(2)cust = 2.0939 . (4.45) 
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4.4 SO(6) x SO(4) model 



Next we evaluate the threshold corrections in the 50(6) x 50(4) string model of 
Ref. ||24jl . The parameters defining this model are also listed in the Appendix. 



This string model realizes the Pati-Salam unification scenario, with gauge group 

S0{6) c x 50(4) ~ SU{A) C x SU(2) L x SU(2) R (4.46) 

and three generations. Here we have only listed that portion of the observable gauge 
group which is of relevance to our analysis, neglecting both the hidden symmetries 
[such as an SU(8)], and various U (l)'s. In the construction of this string model, a ba- 
sis of nine boundary-condition vectors with only periodic and anti-periodic boundary 
conditions is used. In order to achieve the reduction to three massless generations, 
five of these vectors are taken to be the basis vectors of the NAHE set, two others are 
identical to the vectors b 4 and b 5 which appear in the flipped SU (5) model, another 
is similar to the vector 27 of the flipped SU(5) model, and a final vector performs the 
required symmetry breaking from the larger 50(10). Thus, within this construction, 
the 50(6)c gauge symmetry is realized through the L = 1,2,3 fermions listed in 
(|4.1|) , while the 50(4) gauge symmetry is realized through the L = 4, 5 fermions. 



Determining the Charge Polynomials 

Given their similar worldsheet structures, the charge polynomials corresponding 
to the group factors of the 50(6) x 50(4) model are easily obtained from those of 
the first SU(3) x SU(2) x £7(1) model presented above in Sect. 4.2. In particular, the 
three generators of the 50(6) group are simply the generators Qsuts), Qsu(3)i and 
Qu(i) wn i cn appear in (|4.12|) and ( |4.13|) , from which we find that 



Qsom = Qi ■ ( 4 - 47 ) 



Note that all of the non-diagonal terms QlQm with L 7^ M have cancelled in this case. 
Likewise, the generators of the 50(4) symmetry (or equivalently, of the SU{2) L and 

SU(2)r symmetry) are respectively the generators Qsu(2) anc ^ Qu\i) L which appear 
in Q4.15D and (|4.16|) . Their squares are therefore, as before, 

Q 2 SU(2) L = Q\ ~ Q4Q5 

Q 2 su(2) R = Ql + Q4Q5 • (4.48) 



Threshold Contributions: Results of Calculation 

Given these charge insertion polynomials, we see that we have the basis set 
{Qi, Q1Q2, Ql, Q4Q5} of charge polynomials. Our results are as follows: 



37 





boo 


A 




3. 


11.4457 


Q1Q2 


0. 


0. 


g 4 2 


13. 


16.708 


Q4Q5 


-4. 


-3.21395 



(4.49) 



Note that for Q1Q2, the modified partition function B{r) vanishes identically. As 
before, this occurs because every term in B(t) contains a factor of either 0' ^' 
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1/2 



1/2 




, or 0' 
Self- Consistency Checks 



Each of these ©'-functions vanishes. 



We now consider the massless spectrum of this model, and focus on the representa- 
tions of the SO (6) x SO (A) [or S77(4) x SU(2) R x SU(2) L ] observable gauge symmetry. 
The sectors bi 5 produce three generations transforming as (4, 2, 1) © (4, 1, 2) and 
two pairs transforming as (4,1,2) © (4,1,2) under SU{A) x SU{2) x SU{2). The 
Neveu-Schwarz sector produces a pair of Higgs doublets transforming as (1, 2, 2) and 
four sextet fields transforming as (6, 1,1), while the sector S + b 4 + b 5 produces 
an additional pair of Higgs doublets transforming as (1,2,2). The sectors b 14 + a, 
hi + b 2 + b 4 + a, b 2 + b 3 + b 5 + a, and b 4 + b 4 + b 5 + a produce a total of ten 
doublets of SU{2) L and ten doublets of SU{2) R . Finally, the sector 5* + b 2 + b 4 + a 
produces an SU(4) multiplet transforming as (4, 1, 1) © (4, 1, 1). 

The one-loop beta-function coefficients of SU(4:) x SU(2) x SU(2) are given by 



bsu{4) 



lSU(2) L , R 



-12 + - (n A + m) + n 6 



-6 



(4.50) 



where n 4 and n 4 are the total number of massless 4 and 4 representations, where 
n§ is the number of sextet representations, and where n 2i R are the total number of 
doublets of SU(2) LjR respectively. In this model, n 4 + n 4 = 22 (twenty from the 
sectors b 4j 5 and two from the sector S + b 2 + b 4 + a), while n 6 = 4 (all sextets 
from the Neveu-Schwarz sector), and n 2i = 46 and n 2fl = 30. Thus, we find that the 
one-loop beta-function coefficients are 

bsu(4) = 3 , b SU (2) L = 17 , b SU (2) R = 9 . (4.51) 

With the values of 600 given in ( [4.49|) , we see that agreement is again obtained. 

Final Threshold Corrections 
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Given the above results, it is straightforward to calculate the relative threshold 
corrections for each of the Planck-scale group factors in this model, yielding 



&su(2) L - A st/(4) = 8.4763 , A su{2)r - A su{4) = 2.0483 . (4.52) 
We may also consider the threshold correction for the low-energy weak hyper- 



charge. In the Pati-Salam unification scenario [44|, there is an intermediate scale at 



which a further symmetry-breaking occurs, such that 

SU(A) — > SU(3) c xU(1) B -l 
SU(2) R — U(1) R . (4.53) 

Here the U(1)b-l symmetry is generated by the (B — L) generator of the original 
577(4) symmetry, while the remaining U(1)r symmetry is generated by the diagonal 
generator Tj, of the original SU(2)b symmetry. The normalized low-energy weak 
hypercharge is then defined as a linear combination of these two U(l) factors, 

U(l)y = + \Jlu(l)R- (4-54) 

From this relation we deduce that 

2 3 

A U(l)y = g A SU(4) + g A SU(2) R , (4.55) 

yielding 

A u(1)t - A su{4) = - [A SU (2) R - A SU (4)) = 1-229 . (4.56) 
4.5 Non-supersymmetric version of first SU(3) x SU(2) x U(l) model 

In order to assess the role that spacetime supersymmetry might play in affecting 
the overall magnitude of these threshold correction expressions, we have constructed 
a non-supersymmetric version of the first 577(3) x SU(2) x £7(1) model analyzed 
above. This was achieved by altering the GSO projections of that model in such a 
way that spacetime supersymmetry is broken at the Planck scale (e.g., the gravitinos 
are projected out of the spectrum), but no physical tachyons are introduced. Despite 
the absence of spacetime supersymmetry, the spectra of such tachyon-free models 
have a number of interesting properties, such as the appearance of a residual so-called 



"misaligned supersymmetry" [[43], [45] , and the vanishing of certain mass supertraces 
)46[ |. The complete parameters defining this new non-supersymmetric model are listed 
in the Appendix. 

Note that in performing our calculations for this non-supersymmetric model, we 
have continued to use the threshold-correction expressions ( |3.8| ) as originally given 
by Kaplunovsky |T9]. It is claimed in Ref. [BJ| that these expressions hold for all 
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(1, 0) string vacua, regardless of whether or not spacetime supersymmetry is present. 
However, it has recently been shown |5S| that for non-supersymmetric vacua, there 



are also additional terms arising from gravitational interactions such as, e.g., dilaton 
tadpoles.^] Such terms are absent in the supersymmetric case, and do not arise in a 
straightforward "field-theoretic" derivation in which such gravitational interactions 
are not present, or in Kaplunovsky's stringy generalization of such a field-theoretic 
calculation. Nevertheless, they are present in the full string calculation, and are 



found |38j to contribute to the threshold corrections at the same order in a'. In this 
section, however, our goal is merely to construct a toy non-supersymmetric model, 
and to determine how significantly the original expression of Kaplunovsky changes 
when the traces are evaluated over its non-supersymmetric spectrum. Consequently 
we are neglecting those additional terms which would need to be included if we were 
performing a complete analysis of the gauge coupling unification in such a model. 
We hope to perform such a calculation in the future. 

Determining the Charge Polynomials 

Because the boundary conditions of the worldsheet fermions are unaltered relative 
to those of the supersymmetric version of this model, both the gauge group and its 
realization in terms of free fermions remain intact. Therefore the charge polynomials 
for the various group factors are unchanged relative to the supersymmetric case in 
(PH) , fl£T7p , and ([£190 . 



Threshold Contributions: Results of Calculation 

We then repeat the calculation with the same five "basis" charge insertions. Since 
the boundary conditions for all fermions are the same for all sectors as they were in 
the previous spacetime-supersymmetric case, all changes in the results relative to that 
case arise due to the changes in the GSO projection phases. In particular, certain 
particles which had previously been in the spectrum (e.g., various superpartners) 
have been projected out, while new particles which had previously been projected 
out now appear. Our results £1X6 cLS follows: 
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(4.57) 



As in previous cases, the asterisk indicates contributions from unphysical 
tachyons. In fact, the existence of such contributions from unphysical tachyons is 



t We thank E. Kiritsis for discussions on this point. 
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expected in non-supersymmetric string models, since the breaking of spacetime su- 
persymmetry in many cases ensures that the contributions from such unphysical 
bosonic tachyons will not be cancelled by those from any unphysical fermionic tachy- 
onic superpartners. As before, however, the existence of such contributions does not 
lead to any divergence in the corresponding threshold integrals A. We will discuss 
the appearance of these unphysical tachyons in Sect. 6. 

Self- Consistency Checks 

We now turn to the values of 600 • Although the gauge group and charge polynomi- 
als are not altered relative to the supersymmetric version of this model, the breaking 
of spacetime supersymmetry does modify the spectrum. Nevertheless, calculating 
the one-loop beta-function coefficients corresponding to the new non-supersymmetric 
spectrum, we now obtain 

bsu(3) = 1 , b S u(2) = 3 , 

bu(i )o = 17.5 , b u(1)L = 14 , bu(i )Y = 19 . (4.58) 

Comparing this to the appropriate linear combinations of the above values of 600 , we 
see that once again agreement is obtained. 

Final Threshold Corrections 

Using ( |4.14j ), ( f4.17|) , and ( ^4.19|) , we obtain the relative string threshold corrections 
for the SU(3)c, SU(2)l, and normalized U(l)y group factors, 

A u(lh -A sm = 4.934, A u(lh -A sm = 1.9408. (4.59) 

We emphasize once again that these results are to be interpreted as numerical eval- 
uations of the expression given in ( |3.8| ). As discussed above, a complete evaluation 
of the threshold calculations for such a toy non-supersymmetric model would require 



the inclusion of additional terms as well 38 



We see from these results, then, that breaking supersymmetry by a GSO pro- 
jection has no appreciable effect on the magnitude of the expressions (|3 . 8|) . Thus 
the assumption of spacetime supersymmetry in these models is not what sets the 
magnitude of these threshold correction contributions. In Sect. 6 we shall discuss the 
underlying reasons why these string-theoretic threshold corrections are so small. 



5 Confrontation with Experiment 

In the previous section we calculated the threshold corrections Aq within a wide 
class of realistic free-fermionic string models. In this section we will analyze the ef- 
fects of these threshold corrections on the experimentally observed parameters which 
are measured at low energies. In so doing, we shall also consider the other effects 



41 



which might affect the running of the couplings, including non-standard hypercharge 
normalizations, light SUSY thresholds, intermediate gauge structure, and extra non- 
MSSM matter. 

In order to study the effects of the heavy string threshold corrections we have 
calculated in Sect. 4, one method of approach might be to absorb their contribu- 
tions Aq into the logarithms as in (|4.10|) , and thereby obtain a resulting "effective" 



string unification scale. Indeed, this is the course followed in Ref. J39|. One would 
then compare this modified string unification scale with the usual value given for 
^gut; which in turn is obtained by extrapolating the low-energy data under certain 
field-theoretic assumptions. However, we shall instead choose a more general path 
which allows us to take into account a number of different factors which might affect 
the analysis over the whole energy range from the Z scale to the string scale. In 
particular, starting from the string scale, we shall solve the one-loop gauge-coupling 
renormalization group equations (RGE's) for the measured parameters sin 2 6w(Mz) 
and a-^iMz) directly at the Z-scale, and thereby obtain the explicit dependence of 
these parameters on the string threshold corrections. In these RGE's, we shall include 

• the corrections due to second- loop effects. We shall see, in fact, that the two- 
loop contributions are quite sizable, and can in some instances alter our results. 

• the corrections due to Yukawa couplings. 



• the corrections due to scheme conversion. The so-called DR-scheme is used 
in the definition of the string heavy threshold corrections, and in the super- 
symmetric RGE's. However, the MS'-scheme is the one used in extracting the 
low-energy parameters from experiments. One must therefore include explicit 
scheme- conversion terms in the RGE's. 

Furthermore, by explicitly setting up the RGE's in this way, we will also be able to 
include various additional field-theoretic and string-theoretic factors that can affect 
the analysis. These include: 

• the effects due to non-standard values of ki, as can occur in string theory. The 
model-dependent parameter ki is the essentially the normalization of the weak 
hypercharge generator relative to the U(l) generator in a unified theory. In 
SU(5) and 5*0(10) unified models, one has k\ = 5/3. However, in string mod- 
els, the normalization of the U(l) generators is fixed by the requirement that 
the conformal dimension of the massless states be equal to one. The value k\ 
is thus the relative normalization between the properly normalized U(l) gen- 
erators, and the U(l) generator which produces the correct weak hypercharge 
assignment for the Standard Model quarks and leptons. Therefore, depend- 
ing on how these U(l) generators are ultimately realized within a given string 
model, the value of ki may be different from 5/3. 
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• the effects arising from light SUSY thresholds (i.e., the splittings of the sparticle 
mass spectrum). 

• the effects arising from potential additional non-Abelian gauge structure at 
intermediate energy scale. Such additional gauge structure occurs, for example, 
in the Pati-Salam unification scenario, and in the corresponding 5*0(6) x 5*0(4) 
string model. 

• the effects arising from possible additional matter thresholds at intermediate 
energy scales. 

All of these factors contribute, along with the heavy string thresholds, to the running 
of the gauge-couplings, and lead to similar correction terms in the corresponding 
RGE's. Thus, by starting from the string scale and evolving our couplings directly 
down to the Z-scale using such corrected RGE's, we can provide a serious test of 
whether the predictions of string theory are truly in accordance with the experimen- 
tally observed low-energy data. 

Low-Energy Experimental Inputs 

For our subsequent analysis, the input parameters are the tree-level string pre- 
diction [|TJJ 

M s = M string « ^ tring x 5 x 10 17 GeV , (5.1) 



the mass of the Z 47 



M z = 91.161 ±0.031 GeV , (5.2) 
and the electromagnetic coupling at the Z-scale [f47]] 

a' 1 = a e . m .(M z ) -1 = 127.9 ±0.1. (5.3) 

From the RGE's, we then obtain predictions for sin 2 8 w{Mz) and a^(Mz). The 
experimental values we wish to obtain are as follows. For the MS- value of sin 2 ^, 



extracted from low-energy data |48|, we take 



wi 2 e w {M z )\ m = 0.2315 ± 0.001 , (5.4) 



while for the MS- value of as(Mz) we take f49 



ck 3 (M z )Ims = 0.120 ±0.010 . (5.5) 



Renormalization Group Equations 

We seek to compare the above experimentally observed values with the predictions 
from string theory. To do this, we set up our calculation as follows. Recall that in 
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general, string unification implies that each gauge coupling individually has a one- 
loop RGE of the form 

16712 + Mn^% + A? M) (5.6) 



fi'string A* 

where hi are the one-loop beta-function coefficients, and where the /\( total ) represent 
the combined corrections from each of the effects outlined above. From ( |5.6|) we 
wish to obtain expressions for siia 2 9 w(Mz) and a^Mz). To this end, we solve ( |5.6| ) 
for i = 1,2,3 simultaneously in order to eliminate the direct dependence on String 
from the first term on the right sides of (|5.6Q . A small dependence on g str mg remains 
through M string , and in all subsequent numerical calculations we will allow M string to 
vary slightly in order to account for this. In each case we initially assume the MSSM 
spectrum between the Planck scale and the Z scale, and treat all perturbations of this 
scenario through effective correction terms. The expression for a^(Mz) then takes 
the general form 

^\M z )\-ms = aLIsm + Ag + Aff + Ag + 

Afl + AS oop + AW + AW v _, (5.7) 

and likewise for sin 2 0w(Mz)|m^ with corresponding corrections A( sin ). Here A M ssm 
represents the one-loop contributions from the spectrum of the Minimal Supersym- 
metric Standard Model (MSSM) between the unification scale and the Z scale, 
and the remaining A terms respectively correspond to the second-loop corrections, 
the Yukawa- coupling corrections, the corrections from scheme conversion, the heavy 
string thresholds that were calculated in the previous section, possible light SUSY 
threshold corrections, corrections from possible additional intermediate-scale gauge 
structure between the unification and Z scales, and corrections from possible extra 
intermediate-scale matter. Each of these A terms has an algebraic expression in terms 
of well as mo del- specific parameters such as fci, the beta- function coefficients, 

and the appropriate intermediate mass scales. In particular, for sin 2 9yy(Mz) we find: 



A ( sin ) 
^MSSM 



fcl + 1 



l-^(ll-Mln M5 



2tt v ' M 



z 



Ah in) = J-VJ^fo -h )ln M 



M.s 



27T £ h + 1 1 1SP 



sp 



Z 



2tt Y h + 1 V 1 M Mi 

A ( 8 in) = j_ M A 2 -Af 

h ' b ' 21^ + 1 2 v 7 

where M s = M string is the string unification scale, a = a e . m .(M^), M sp are the spar- 
ticle masses, Mj are the intermediate gauge mass scales, and A 3 2 Y are respectively 
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the heavy string threshold corrections for the strong, electroweak, and properly nor- 
malized hypercharge gauge group factors. Likewise, for (Mz), we find: 



A («) 
^MSSM 



A 



(a) 
Is. 



" + ^-(-3*i -15) In ^ 
a lit Mz 

1 



A 



(«) 



A 



(a) 
h.s. 




, bo — bi 

1 + h 1 sp sp 



1 



1 + h 



K - b 3 



M z 
M s 



In 



Mi 



4tt 1 + h 



[h(Ay - A 3 ) + (A 2 - A 3 )] 



(5.9) 



Note that in the above solution of the one-loop RGE's for sin 2 9w(Mz) and 
«3 (Mz), only the differences of the RGE's ( |5.6|) for the separate gauge couplings 
are used. Therefore, the expressions for siia 2 9 w(Mz) and as(M z ) depend only on 
the differences of the heavy string threshold corrections A 3 2 Y of the different group 
factors, and not on their absolute values. Consequently, the group-independent ad- 
ditive factors in the full expression for the one-loop string threshold corrections do 
not affect the predictions for the experimentally measured parameters sin 2 9w(Mz) 
and a^(M z ). This is important, for these extra factors are neglected in the original 
definition for Aq in (^H). 

Finally, before we can proceed, we must estimate the second-loop and Yukawa- 
coupling corrections. As discussed above, these are calculated assuming the MSSM 
spectrum between the string unification and Z scales. To estimate the size of the 
second-loop corrections, we run the one- and two-loop RGE's for the gauge couplings 
and take the difference. Likewise, to estimate the Yukawa- coupling corrections, we 
evolve the two-loop RGE's for the gauge couplings coupled with the one-loop RGE's 
for the heaviest-generation Yukawa couplings, assuming X t w 1 and X b = \ T rs 1/8 
at the string unification scale. We then subtract the two-loop non-coupled result. 
As indicated above, these differences are then each averaged for different values of 



M 



string 



(in order to account for various values of String )• Numerically, this yields the 



following results for sin 9w(Mz)\jys : 



A ( sin ) 
AA 2-loop 

A ( sin ) 
^Yuk. 
^ (sin) 



while for a 3 1 (Mz)\jy s - we find: 



a(o) 

AA 2-loop 
^Yuk. 

a (q) 



0.00678 

-0.00091 
0.00026 , 



0.80757 

-0.16728 
0.0596831 . 



(5.10) 



(5.11) 
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In the next subsections, we will discuss and evaluate the various remaining con- 
tributions to the experimentally measured parameters sin 2 6 ' W (M Z ) and a 3 (M z ). 



5.1 MSSM spectrum between M z and M string 

In this section we discuss the predictions for sin 2 6\y(Mz) and a 3 {Mz) assum- 
ing that the particle spectrum between the string unification scale and the Z-scale 
is purely that of the Minimal Supersymmetric Standard Model (MSSM). The con- 
tributions of such a spectrum are those given in the first lines of ( |5.8|) and ( |5.9| ) 
respectively We also include the two-loop and Yukawa-coupling corrections to the 
one-loop RGE's, and the corrections due to the conversion from the DR-scheme to 
the MS'-scheme. These are the corrections which appear in ( |5.10|) and ( 5.11|) . 



Using the RGE's with these two terms as discussed in Sect. 5.0, we find that 
the prediction for sin 2 6 ] w(Mz) deviates from the experimentally measured value by 
approximately 5%, and that the prediction for a 3 (Mz) deviates from the experi- 
mentally measured value by approximately 100%. The discrepancy with a 3 (M z ) is, 
of course, expected, since it is this discrepancy which is ultimately the root of the 
factor of twenty that separates the string scale M string from the MSSM unification 
scale M M ssm- Indeed, in this way we can even see why the value of a 3 (M z ) should 
be off by 100%. Suppose, for example, that we were to run the RGE's not from 
^string ~ 20Mmssm to Mz, but rather from Mmssm to M^/20. As expected from the 
running of the QCD coupling below the Z scale, we would then find that a 3 {Mzl '20) 
deviates from a 3 (M z ) by approximately 100%. However, due to the additive nature 
of the logarithms involved, this effect on 0:3 should be identical to that achieved by 
running instead from M str i n g to Mz- Indeed, in the string unification scenario, we 
have simply shifted this discrepancy from the low scale near Mz to the high scale 
near M string . Thus, we see that the simplest string unification scenario assuming the 
MSSM spectrum below the string scale should indeed predict a value of a 3 [Mz) that 
is off by 100%. 

In Figs. H and ||, we have plotted sin 2 9w(Mz) and a 3 {Mz) versus M str i ng in the 
range 1 x 10 16 GeV < M str i ng < 7 x 10 17 GeV, both with and without the two-loop 
corrections. In the analysis for these figures, rather than using the averaged second- 
loop corrections from ( 5.1 U| ) and ( p. 11 ), we have run the full two- loop RGE's over 



the entire range of M str i n g- In Fig. ||, we have plotted the two low-energy variables 
sin 2 9w(Mz) and a 3 (Mz) versus each other as M string is varied. From these figures we 
see that for 3 x 10 17 GeV < M string < 7x 10 17 GeV, both sin 2 6 W {M Z ) and a 3 (M z ) are 
approximately ~ 0.22. Thus, as expected, assuming the MSSM spectrum between 
the Z scale and the string unification scale results in significant disagreement with 
the low-energy data. 
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Figure 1: Dependence on the string scale: sin 2 9w(Mz) vs. M f 
one-loop and two-loop running are plotted. 



string 



Results for both 



5.2 Non-standard hypercharge normalizations 

We can also examine the extent to which this MSSM string unification scenario 
depends on the value of the hypercharge normalization k\. In the 577(5) or 50(10) 
unification schemes, the value of k\ is 5/3. However, in string unification, this as- 
sumption is not necessary, as other embeddings of the weak hypercharge into the 
four-dimensional string gauge group are possible. This possibility is also discussed 
in Ref. JL8|. In Figs. f| and |5j, we plot the predictions for sin 2 6 w (Mz) and a 3 (M z ) 
versus k x in the range 1 < k x < 2. We take 3 x 10 17 GeV < M string <7x 10 17 GeV. 
From these figures we find that the one-loop string prediction can be in agreement 
with the measured values for a value of k\ w 1.4. 

Unfortunately, none of the realistic models we have been analyzing has a value 
of k\ of this size, and indeed all of the models have k\ > 5/3. Similar results have 
also been found, for example, in the case of semi-realistic orbifold models (see, e.g., 
Ref. ||50|| ). Moreover, in string theory, there are general reasons to expect that k\ > 
5/3 in realistic models. The reason is as follows.f] As discussed at the end of Sect. 4.3, 
the weak hypercharge is a combination of simple worldsheet currents which are each 
normalized to one. To produce the correct conformal dimension for the massless 
states, every £7(1) generator (each of which is ultimately a combination of simple 
worldsheet currents) must be normalized to one. By contrast, the U(l) generator 
that produces the correct weak hypercharges for the Standard Model particles is not 

* We thank J. March- Russell for discussions on this point [ jsi) . 
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Figure 2: Dependence on the string scale: a^Mz) vs. M string . Results for both 
one-loop and two-loop running are plotted. 

normalized to one, and k\ is essentially the normalization coefficient of the properly 
normalized weak hypercharge generator. Since the weak hypercharge generator is a 
combination of simple worldsheet currents, we can determine the minimal number 
of simple worldsheet currents that must be used in order to generate the correct 
hypercharges for all the quark and lepton families and to satisfy the various anomaly 
cancellation constraints. We then find that the minimal value of k\ is essentially 
5/3. It is not likely, therefore, that this effect can explain the discrepancy between 
the low-energy data and the suppositions of string unification. Thus, in the analysis 
below, we shall take k± — 5/3 unless otherwise stated. 

5.3 Light SUSY thresholds 

We now study the effect of the light SUSY thresholds on the predictions for 
sin 2 Q w {Mz) and a 3 (M z ). These effects correspond to the second lines of (|5.8| ) and 
( [5.9|) respectively. 

Our purpose here is not a detailed quantitative analysis of the sparticle threshold 
corrections, but rather a qualitative examination of whether the light SUSY thresh- 
olds are capable of removing the discrepancy, found in Sect. 5.1, between the predicted 
and the experimentally observed values of sin 2 9w(M z ) and a 3 (M z ). Consequently, 
to parametrize the sparticle thresholds, we shall consider the spectrum of the MSSM 
near the Z scale, and neglect the contributions from the Yukawa couplings and elec- 
troweak VEV's. This implies that all of the supersymmetric scalar mass matrices, 
including that of the stop quark, will be diagonal, and that the D-term contributions 
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Figure 3: Dependence on the string scale: a 3 (Mz) vs. sin 2 9w(Mz), as M str mg is 
varied. Results for both one-loop and two-loop running are plotted. 



to the sparticle masses are neglected. With these assumptions, the sparticle masses 
can be obtained from the one-loop RGE's for the soft SUSY-breaking terms. In terms 
of the usual soft SUSY-breaking parameters m , m^, and /i, the gluino, wino, and 
higgsino masses are then given by M 



m? 



m 1/2 



a 2 {m § 
rriyj = — mi/2 



a 3 a 3 
Likewise, the scalar sparticle masses are given by 



mr =fj, 



m l = m l + Cp m l/2 



(5.12) 



(5.13) 



where the coefficients cp for the different sparticles are given in terms of their hyper- 
charges Yp by 

c 3 (m p ) + c 2 (mp) + Y~c y (m p ) 



p 



with 



c 3 (m p ) 
c 2 {mp) 



c y {mp) 



9 

3 
2 

2h r 
TT 



■p ^y\""pj 
1 - (l + 3X)~ 2 

1 - (1 -x)- 2 
i - (i - nx/hy 2 



(5.14) 



and 



X 



1 i -^string 

~Z ^string m 

2tt ma 



(5.15) 
(5.16) 
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Figure 4: The effects of hypercharge normalization: sin 2 6 'w(Mz) vs. k\. The width 
of the curve reflects the variations in M string . 

The corresponding corrections to sin 2 Ow(Mz) and a^Mz) due to these light SUSY 
thresholds are then given by 

(sin) _ J_ ha v r (sin) , run 

and 

A[ a J = —y^n R CP\n^ (5.18) 

where the summations are over the particle representations R with masses uir. Here 
hr is the number of times a particular generic representation appears in the complete 
spectrum {e.g., three for three generations), and the C^-coefficients are given in terms 
of the one-loop beta-function coefficients b 32 f(R) of the representation R as 

4 Sill) = by(R)-b 2 (R) 

4 a) = ^by(R) + ^b 2 (R) - b 3 (R) . (5.19) 
Note that b Y (R) = b Y {R)/ki. For each representation R, these beta-function coeffi- 
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Figure 5: The effects of hypercharge normalization: a^{Mz) vs. k\. The higher curves 
correspond to higher values of M string . 



cients, degeneracies ur, and C^-coefficients are as follows: 
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Assuming universal boundary conditions for the soft SUSY-breaking terms {i.e., 
assuming that m 1( / 2 and mo are universal), we have analyzed the possible light 
SUSY threshold contributions for a wide range of points in the parameter space 
{mo,mi/2,m(i,mJ. In general we find that the light SUSY thresholds are small, 
and that the predictions for sin 2 6 ' W (M Z ) and a 3 (M z ) continue to disagree with the 
measured values. In Fig. | we display sin 2 6{Mz)jis versus a 3(^z)is f° r a sampling 
of points in the SUSY-breaking parameter space. Within this parameter space, each 
free parameter X is sampled in the interval (X min , X max ) with spacing AX between 
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consecutive points, as follows: 



Parameter X 



m (GeV) 
mi/2 (GeV) 
m h (GeV) 
m~ h (GeV) 
(GeV) 



string 



a. 



-1 
string 



V m ; n 
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100 


500 


200 
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7 x 10 17 


5 x 10 16 


20 


32 


4 



(5.21) 



We have taken the top-quark mass m t = 175 GeV. 
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Figure 6: Scatter plot for {sin 2 6 W (Mz), a 3 (M z )} for various values of {m ,mi/2, 
m h , m~ h , M string }, assuming the MSSM spectrum. 

It is important to note that our assumption of universal boundary conditions 
for the soft SUSY-breaking terms is not necessary. In fact, in string models, such 



boundary conditions are generally expected to be non- universal |52] . We can then ask 
whether non-universal boundary conditions can bring the string-scale unification into 
agreement with the experimentally measured couplings at low energy. By examining 



the different terms in ( |5.17| ) and ( |5.18 ), we can see that for this purpose, the value 



of the light SUSY threshold corrections must be positive. Now, the signs of the Cr- 
coeflicients of the different terms in these two equations are determined by the one- 
loop beta-function coefficients of the different representations R. Assuming that all 
the sparticle masses are larger than Mz, we can maximize the light SUSY corrections 
by setting the masses of the negative contributions equal to the Z mass. For example, 
in ( |5.18|) , we see from the signs of the coefficients that the wino, top, sleptons, 
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higgs and higgsino all give positive contributions. Assuming a common mass M 
for all of these but the top, and summing over the resulting coefficients, we obtain 
A£? = 15/(16tt) In (M/M z ). Thus, with a <f redicted ) „ .2 and a (~ red ) ~ .14, even 
the best-case scenario requires M rs 100 TeV. For sin 2 6w, on the other hand, only u r , 
d r , and e r contribute with positive coefficients. Assuming that all the other thresholds 
are degenerate at Mz, and that u r , d r , and e r have a common mass M, we find that 
agreement with the measured value of sin 2 9w requires M ~ 10 TeV. Alternatively, 
some of the sparticle thresholds may be below Mz- Assuming degenerate thresholds 
M = 45 GeV for the sparticles making negative contributions to sin 2 9w, and with 
m t = 175 GeV, we find that the contribution for Aj^ from these thresholds is of 
the order of 0.001, which still requires M ~ 10 TeV. Thus, even with these extreme 
"best-case" scenarios, light SUSY thresholds cannot by themselves bring the string- 
unification predictions into agreement with the experimentally observed values. 

5.4 Heavy string thresholds 

We now analyze the effects of the threshold corrections due to the infinite towers 
of heavy string states. These corrections were the focus of our calculations in Sect. 4, 
and contribute to the running of the gauge couplings as indicated in the fourth lines 
of (U) and (p|). 

The calculation of the heavy string threshold corrections in Sect. 4 was performed 
for normalized gauge group factors Cr, with ki = 1. However, as indicated in (|3.3|) , 
the full string threshold corrections are in general of the form Aj = Aj + kiY, where 
Y is independent of the gauge group. In Sect. 4 we calculated the gauge-dependent 
terms Aj for the different gauge group factors. Therefore, by taking the differences of 
the threshold corrections for two different gauge groups, the group-independent term 
cancels. Moreover, as we see from (|5.8|) and ( |5.9|) , it turns out that the solutions 
of the RGE equations for sin 2 8 ' W (M Z ) and a 3 (M z ) depend only on the differences 
Aj — Aj for the different group factors Gj and Gj in a given string model. There- 
fore, the low-energy predictions for sin 2 6w(Mz) and a^Mz) do not depend on the 
group-independent factor Y . One might worry, of course, that the presence of a 
non-standard normalization for k\ in certain string models might render this claim 
incorrect. However, the calculation of the threshold corrections is always done with 
respect to the properly normalized U(l) generators. Similarly, only the threshold cor- 
rections for the properly normalized U(l) generators appear in the RGE equations. 
This is similar to the usual practice in grand unified theories in which the coupling 
of the properly normalized U{1) generator evolves below the GUT scale, and the 
matching with the weak hypercharge coupling is done at the Z scale. Therefore the 
value of the gauge-independent term Y continues to be irrelevant for our analysis. 

We now discuss the effects that these string threshold corrections will have on the 
low-energy parameters. Because the gauge structure and threshold corrections are 
highly model-dependent, we shall have to consider each realistic free-fermion model 
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in turn. 

Considering first the string model analyzed in Sect. 4.2, we see that the weak 
hypercharge has the standard 5*0(10) embedding. Therefore, we can directly insert 
the relative values for Asu(3), Asu(2) and Ajj(i) t that we found in (|4.23|) into (|5.8| ) 
and ( |5.9| ). For this model, we consequently obtain 

a£J° = -0.0006 and A n Q s } = -0.3536 . (5.22) 

Thus, the effect of the string threshold corrections in this model is to slightly decrease 
the values of sin 2 6w{M z ) and a^ 1 (M z ). Alternatively, of course, this can be regarded 
as effectively increasing the string unification scale. This can be seen by absorbing 
the string threshold corrections into ln(M string /M^) in ( |5.8| ) and (|5.9|). From 
we see that the corrected string unification scale can be written as 



, ^(corrected) _ , , (uncorrected) 

m string ~~ Jw string ex P 



fci(Ax- A 3 ) + (A 2 - A 3 
2(3/c! + 15) 



(5.23) 



Numerically, this yields 

M^Zt^ ~ 6 - 72 x 1017 GeV ( 5 - 24 ) 

where we take Mg™™ rrectcd ' ) « 5 x 10 17 GeV. Thus, the effect of the string threshold 
corrections in this model is to enhance the disagreement with the experimentally 
observed values. 

We next analyze the model of Sect. 4.3. As we saw in Sect. 4, in this model the 
analysis is complicated due to the presence of the enhanced symmetry. The group 
factors relevant for the analysis are SU(3)c X SU(2)l X SU(2) CVLSt x U(1)y>- Above 
the S77(2) cust x U(1)y> breaking scale, we have to consider the running of these four 
group factors. The S77(2) cust x U(1)y' symmetry can be broken down to E/(l)y by, 
for example, the VEV of the right-handed neutrino, and below this breaking scale 
we have the standard one-loop RGE's for SU(3)c X SU{2) L x U(1) Y . The relation 
between the weak hypercharge coupling and the SU(2) CVLSt x C/(l)y couplings is then 

1 1 17 , 

(5.25) 



«y 4a cust 12 a' 



Y 



where a.y t is the coupling of the properly normalized U{1) Y ,. The RGE for the U(l)y 
coupling then takes the form 

1 3 1 3b Y Mj I 3.1 17, M s 

m —-(-02c + T7:0v') m - 



"string 5«i(/i) 10vr n 2tt5 v 4 12 Y ' Mj 

Let us assume for the moment that the VEV of the right-handed neutrino is of the 
order of Ms- This implies that Mi = Ms- Therefore, below the string scale we have 
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the three gauge couplings of the Standard Model, and the threshold corrections for 
the properly normalized weak hypercharge are given by ([4.41 ) . Inserting the relative 



values of A 3 , A 2 , and Ay from ( f4.41p into ( |5.8j ) and ( |5.9| ), we then obtain 



A[ si s n) = -0.0012 and A£J = -0.335 . (5.27) 

As in the case of the previous model, we see that the effect of the string threshold 
corrections is to decrease the values of sin 2 9w(M z ) and «3 1 (M^), or equivalently 
to increase the effective string unification scale. Using ( p.23[ ) with Mj t ™° rrected ' « 
5 x 10 17 GeV, we now obtain 

^string eCted) « 6.62 x 10 17 GeV . (5.28) 

We shall shortly discuss the modifications that arise if Mj < Ms- 

Next, we analyze the effect of the threshold corrections in the 5*0(6) x SO (4) string 
model which realizes Pati-Salam unification scenario. In this string model, there is 
an intermediate energy scale at which the SU{4) x SU(2)l x SU(2)r symmetry is 
broken to SU(3)c X SU{2) L x U(l)y- The low-energy parameters are consequently 
affected by the presence of this intermediate energy scale. The weak hypercharge Y 
is a linear combination of B — L and T^ R : 




Y = n TB - L + w 3 ' ■ (5 - 29) 

The low-energy parameters are then obtained from the one-loop RGE's via 

1 1 b 3 Mj 6 4 . M s A 4 
H In — — I m — — H 



a 3 (M z ) a string 2tt M z 2tt M/ 4tt 

1 1 b 2r , M 7 b 2l , M s A 2r 

+ In — + In — + 



a 2 (M z ) a string 2vr M z 2vr M 7 4tt 

1 1 &! Mj 1 2 3 M s 

+ — ln^-^ + —(-64 + -0 2 „)ln^- 



cn(M z ) a string 2vr M z 2tt v 5 5 R> Mj 

+ ^(^A 4 + ^A 2 J (5.30) 

Assuming again that M/ « M5, the threshold corrections for the three properly 
normalized group factors are A 4 , A 2i and A y = (2/5)A 4 + (3/5)A 2fl . Inserting the 
values from ( |4.52j ) and (|4.56|) into ( p.8| ) and (|5.9|), we find 



A[ sl s n) = -0.0028 and a£ = -0.3141 (5.31) 

Just as in previous models, the effect of the string threshold corrections is therefore 
to increase the effective string unification scale. Inserting the values of A 4 , A 2l and 



Ay into ( |5.23| ), we obtain 

Mj t c °" ected) « 6.5 x 10 17 GeV . (5.32) 



55 



Thus, in all of the examples that we have explored, we find that the effect of 
the string threshold corrections is to increase the effective string unification scale. 
Therefore, the effect of the string threshold corrections in these examples is always 
to enhance the disagreement with the low-energy observables. 

It is an important observation that the sizes of the heavy string threshold correc- 
tions are very small in all of these realistic string models, and thus do not greatly 
affect (either positively or negatively) the magnitude of the string unification scale. 
Although such threshold corrections receive contributions from infinite towers of mas- 
sive string states, we have in fact been able to provide a general model-independent 
argument which explains why these corrections are naturally suppressed in string 
theory (except of course for large moduli). This will be discussed in Sect. 6. Thus, 
we conclude that these threshold corrections cannot by themselves resolve the exper- 
imental discrepancy. 

5.5 Intermediate gauge structure 

In the last two models that we studied at the end of Sect. 5.4, there exists an 
intermediate energy scale Mj at which an extended gauge group is broken to the 
gauge group of the Standard Model. In each case, we assumed that the extended 
symmetry is broken at the string scale, and that therefore Mi = Ms- An obvious 
issue, then, is to determine the effect of the extended symmetry if it is broken at an 
intermediate scale Mi < Ms- In particular, might the breaking at the intermediate 
scale bring the string-scale predictions into agreement with experiment? 

We claim that in the above examples, this cannot happen. Our reasoning is as 
follows. Since the larger discrepancy with experiment is for a^, it is sufficient to focus 
on this observable. Now, in the first model we considered, the RGE for a 3 does not 
depend on the intermediate scale breaking at all. Thus the disagreement with the 
experimentally observed value persists for all Mi. Of course, this argument is some- 
what deficient due to the dependence of the RGE for on the group-independent 
contribution Y to the heavy string threshold corrections. However, unless (and con- 
trary to expectation |k|) this group-independent term Y is very large, the argument 
will hold. A more careful argument is provided by solving the RGE's, in the presence 
if intermediate scale Mi, for sin 2 9 w(Mz) and a^(Mz)- The one-loop RGE for a.\ 
takes the form 

JL = JL_ _ L^Mf _ ± Af + J. 1 {bY _ A + ln^ (5.33) 

«i k x a Y 271k! M z 4vr y 27rV v 4 12 M I y 1 

where the first three terms in ( |5.33|) are identical to the RGE without an intermediate 
scale, and the last term incorporates the effect of the intermediate scale. Thus, in 
the presence of the intermediate scale the solution of the RGE's for sin 2 9 w(Mz) in 
(15.81) has an additional term 



1 a 
+ — 



2vr fci 



by ~ ( T + T2 b YI ) 



mf" (5.34) 
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Since sin 2 9w appears divided by a = a e . m . in the solution of the RGE's for 0:3 , the 
effect of the intermediate scale on (Mz) is that in ( |5.34| ) divided by a. Now, for 
a spectrum consisting of three generations and two Higgs doublets, the coefficient 
of the logarithm in ( 5.34j) is —301/48. The crucial point here is the sign of this 
coefficient, as it decreases both sin 2 9 'w(Mz) and a 3 ~ 1 (M^). Consequently, a^{Mz) 
increases, and the disagreement with experiment is enhanced. 

Similar considerations also apply to the 50(6) x SO (A) model. From Q5.30D , we 
see that the effects of the intermediate scale on sin 2 9 W [M Z ) and a 3 (M z ) amounts 
respectively to the additive terms 

4 ha M s ( . 

- + AMTTki^M; (5 - 35) 

and 

1 - M s 




= ... - — In— - (5.36) 

in ([5.8D and (|5.9|) . We assume that the spectrum below Mj is that of the MSSM, 
and that above Mj it consists of three 16's of 5*0(10) that produce the three genera- 
tions, one (1,2,2) representation that produces the light Higgs, and (4, 1, 2) + (4, 1, 2) 
representations that are used to break the SU(4) x SU(2)r symmetry at Mi. The 
signs of the coefficients in (|5.35|) and (|5.36|) shows that the effect of having Mj < M s 
is to increase sin 2 9 w(Mz) and a^{Mz)- Note that in arriving at this conclusion, 
we have used only the differences of heavy string threshold corrections, which are 
unambiguous. Thus, the effect of having an extended gauge structure broken at 
an intermediate energy scale is to enhance the disagreement with the experimental 
results for a 3 (Mz). 

Likewise, in the case of the flipped SU(5) model, the effect of the intermediate 
scale can be incorporated via the following additional term in 



1 32 ha M s 
~ 27 IT k^+l M~j ' (5 ' 37) 

Here we have assumed the spectrum below Mj to be that of the MSSM, and above 
Mj to consist of three 16 representations of 50(10), one 5 and 5 of SU(5) that 
produces the light Higgs doublets, and one 10 and 10 of SU(5) that is used to 
break the SU(5) x 17(1) symmetry to SU(3) x SU{2) xU{l). Thus, the effect of the 
extended gauge structure in this model is to reduce sin 2 9 W (M Z ), and to enhance the 
disagreement with the experimentally observed value. 

It is remarkable that in all of these realistic string models, the availability of an 
intermediate scale Mj does not help in removing the discrepancy between string- 
scale unification and low-energy data. A priori, one might have expected that the 
presence of such an extra degree of freedom would have enabled agreement to be 
reached. However, we now see that within the context of the realistic string models, 
this extra degree of freedom only worsens the agreement. Thus, we can effectively 
rule out intermediate gauge structure as a potential explanation for the discrepancies. 
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5.6 Intermediate matter thresholds 



We now turn to the effects induced by additional matter below the string scale 
17, 54| . As we shall see, such matter appears naturally in the string models we 



have examined. It is evident from our analysis up to this point that such intermediate 
matter now appears to be the only way by which the disagreement with experiment 
can be possibly be resolved. Fortunately, as we shall find, in certain models exactly 
the required matter appears, in just the right representations and with just the right 
non-standard hypercharge assignments. Hence, in these models, the disagreement 
with experiment can be resolved. 

In the analysis up to this point, we have assumed that the matter spectrum 
below the string scale is simply that of the Minimal Supersymmetric Standard Model, 
i.e., that it consists of exactly three chiral generations with two Higgs doublets. 
However, in the context of the realistic string models, this assumption is ad hoc. In 
fact, in all of the string models constructed to date, additional color triplets and 
electroweak doublets beyond the MSSM appear in the massless spectrum, in vector- 
like representations. The number of such additional color triplets and electroweak 
doublets is of course highly model-dependent, as are their mass scales. However, the 
assumption that the massless spectrum below the string scale is that of the MSSM 
is, in general, not justified. 

Mass terms for these extra states beyond the MSSM may arise from cubic or 
higher-order non-renormalizable terms in the superpotential. In the models studied 
to date, the mass scale of the additional color triplets and electroweak doublets is not 
the string scale. In general, the masses of the extra states are suppressed relative to 
the string scale, the reason being that the mass terms arising from non-renormalizable 
terms are suppressed relative to the cubic- level mass terms. For example, in Ref . [^5|] . 
the mass scale of an additional pair of color triplets was estimated to be of the order 
of 10 11 GeV. 

The contributions of such additional color triplets and electroweak doublets to 
the low-energy parameters are given in the third lines of ( |5.8| ) and ( |5.9| ). Provided 
that these additional states exist at the appropriate scales, we shall find that the 
presence of this additional matter results in agreement between the hypothesis of 
direct unification at the string scale, and the values of the low-energy parameters 
sin 2 Q W {M Z ) and a 3 (M^). Of course, there exist a large number of possible scenarios 
for the mass scales of the extra states which will allow direct unification at the string 
scale, and a classification of all these possibilities is beyond the scope of this paper. 
Indeed, such an analysis would require a detailed investigation of each individual 
string model, its full spectrum, and its renormalizable and non-renormalizable su- 
perpotentials. This is currently being investigated [Q. However, the essential point 
that cannot be overemphasized is that these string models, in general, produce just 
the sort of spectrum with additional matter as required to have the possibility of 
direct string unification. 
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In the realistic free-fermionic models we have studied, the three generations are 
obtained from the sectors b 1; b 2 , and b 3 . The electroweak doublets that couple to 
the states from these sectors, and which therefore may correspond to the MSSM 
Higgs doublets, are obtained from the Neveu-Schwarz sector and a sector which is a 
combination of the vectors {bi, b2, b3, a, /?}. For example, in the models of Sects. 4.2 
and 4.3, the electroweak doublets are obtained from the Neveu-Schwarz sector and 
the bi + b 2 + a + (3 sector. The Neveu-Schwarz sector and the b x + b 2 + a + (3 sector 
may, in general, also produce color triplets. Color triplets from these sectors couple to 
the massless fermions from the sectors b 1; b 2 , and b 3 . Their mass scale is therefore 
restricted from proton decay through Higgsino exchange. However, for these two 
sectors there exist a superstring doublet-triplet splitting mechanism through which 
the color triplets are removed from the spectrum via GSO projections M. 

In addition to the massless spectrum from the above sectors, the free-fermionic 
models may contain additional color triplets and electroweak doublets form several 
other sectors. In the SO(6) x SO(4) and the flipped SU(5) x U{1) models, additional 
16 and 16 representations are obtained from the additional vectors beyond the NAHE 
set that are used to reduce the number of generations to three. Typically, there are 
one or two of such pairs. These produce the (3, 2) 1( / 6 , (3, 1)1/3, (3, l)-2/3, and (1, 2)_i/ 2 
representations of SU(3) x SU{2) x £/(l)y. These states have the usual one-loop 
beta-function coefficients. In the standard-like superstring models, additional color 
triplets and electroweak doublets may also appear from sectors which arise from 
combinations of the vectors {bi, b 2 , b 3 } + {a, (3, 7}. These states, in general, do not 
have the standard weak hypercharge assignments, and therefore their one-loop beta- 
function coefficients will be different from those of the MSSM representations. For 



example, in the model of Ref. [17|, there exist color triplets {-D3, D 3 } and electroweak 



doublets {£, £} with the following beta-function coefficients: 

( \ 

(5.38) 
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To estimate the effect of the intermediate thresholds induced by such additional 
matter, we shall examine in detail the spectrum of the model of Sect. 4.3. In order to 
be as general as possible, we will start from (|5.8|) and (|5.9|) and from the experimental 
constraints on the low-energy observables, and derive the corresponding constraints 
on these intermediate matter thresholds. As our low-energy experimental constraints, 
we shall impose 

0.230 < sin 2 9 W (M Z ) < 0.233 

0.110 < a 3 (M z ) < 0.135 . (5.39) 

As before, we set M string = 5 x 10 17 GeV, and take k x = 5/3 and a e . m .(M z ) = 1/127.9. 
We set all of the sparticle masses to be degenerate at M z . These values, along with the 
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corrections from two-loop contributions, Yukawa couplings, and scheme- conversion, 
as well as the values for Aj given in ( |4.41|) , are then inserted into the expressions for 
sin 2 Q w {Mz) and a 3 (M z ) given in ( |5.8| ) and ( p.9[ ) respectively. Using our experimental 
limits ( |5.39|) , we then obtain the constraints 



M<? 

10.29 < Zi(b 2i -b h ) ln^ < 14.14 



39.22 < 



hb u + b 2i - (1 + fc0& 3 



In ^ < 67.40 . 

Mi 



We can also combine these two equations to obtain 



18.57 < £(& 3i - 6lj )ln|| < 30.55 



(5.40) 



(5.41) 



We can now use (|5.40|) and (|5.41|) to examine the possible general scenarios for 
intermediate matter states. Our first observation is that in order to accommodate 
the low-energy parameters, both intermediate color triplets and electroweak doublets 
are needed. Let us now focus on the representations that are available in the string 
model of Refs. [|T7], |23"| . In this model, there are three color triplets in vector-like rep- 
resentations, two with the one- loop beta-function coefficients of (|5.44|) , and one with 
the one-loop beta-function coefficients of ( |5.38| ). There are five pairs of electroweak 
doublets from the Neveu-Schwarz and b x + b 2 + a + (3 sectors with the quantum 
numbers 

( b SU(3)\ ( \ 

6*1,(2) = 1/2 , (5.42) 



and three pairs of electroweak doublets with the quantum numbers of (|5.38|) . Let us 



suppose that all of the intermediate thresholds are from states that fit into the 5 and 
5 representations of SU(5). For one such pair, we obtain from ( |5.40| ) and ( |5.41| ) the 
constraints 



30.95 < In 



27.725 < In 
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< 



(5.43) 



where the upper limits are not important. It is clear that these two equations cannot 
be satisfied for any values of M 2 and M 3 , or for any number of doublets and triplets. 
We therefore conclude that in any string models containing only those states that 
fit into 5 and 5 representations of SU(5), intermediate matter thresholds cannot 
account for the disagreement with the experimentally observed values of sin 2 0w(Mz) 
and a^{Mz)- It is an important conclusion that some string models can actually be 
ruled out in this way. 
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By contrast, in the model of Sect. 4.3, the massless spectrum contains not only 
states that fit into the 5 and 5 representations of £77(5), but also color triplets and 
electroweak doublets with exotic weak hypercharge assignments. These therefore 
cannot be fit into £0(10) representations. In this particular model, there are two 
pairs of color triplets {D\, D\, 7J 2 , D2} with one-loop beta-function coefficients 

(bsu(3)\ fl/2\ 
bsu(2) = (5.44) 

V b U(l) J DU D U D 2 ,D 2 V 1 / 5 / 

in addition to one pair of color triplets with the quantum numbers of ( |5.38| ), and 
three pairs of electroweak doublets with the quantum numbers of Q5.44j ). We shall 
set the masses of these three doublet pairs to be degenerate at one scale, M 2 , and 
examine possible scenarios for the masses M 3 of the color triplets. With one light 
color triplet, for example {D\,Di}, we obtain the limits 

experimental limit < M 3 < 18141 GeV . (5.45) 

Setting M 3 at the upper limit of ( |5.45| ), we find 

7.2 x 10 13 GeV < M 2 < 2.6 x 10 14 GeV , (5.46) 
while for a lower limit of M 3 ~ 500 GeV we obtain 

3.6 x 10 5 GeV < M 2 < 1.7 x 10 6 GeV . (5.47) 
By contrast, with two triplet pairs degenerate at one mass scale M 3 , we instead find 

4.3 x 10 6 GeV < M 3 < 9.5 x 10 10 GeV , (5.48) 
so that taking the upper limit for M3 yields 

7.2 x 10 13 GeV < M 2 < 2.6 x 10 14 GeV (5.49) 
while the lower limit on M 3 yields 

5 x 10 12 GeV < M 2 < 1.8 x 10 13 GeV . (5.50) 
Finally, with all three color triplet pairs degenerate at the scale M 3 , we find 

2.4 x 10 11 GeV < M 3 < 7.2 x 10 13 GeV (5.51) 

for which the upper and lower limits respectively yield 

3.7 x 10 14 GeV < (M 2 ) uppcr < 1.1 x 10 15 GeV 

5.7 x 10 13 GeV < (M 2 ) lowcr < 2 x 10 14 GeV . (5.52) 
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Clearly, many viable scenarios exist, and the above examples are not exhaustive. 
In Fig. we plot sin 2 9 W (M Z ) versus a 3 (M z ), with M 3 = 2 x 10 12 GeV and M 2 = 
2 x 10 14 GeV. We have included in this analysis the two-loop, Yukawa, and scheme- 
conversion corrections given in (|5.10|) and ( 5.11 ), as well as the heavy string threshold 



corrections and the light-SUSY corrections due to the splitting of the sparticle spec- 
trum. The unification scale is varied between 3 x 10 17 GeV < M str i ng < 7 x 10 17 GeV, 
and 0.03 < a str m g < 0.05. It is evident from the figure that in this model, the low- 
energy experimental parameters can indeed be accommodated, provided that the 
thresholds from this extra matter exist at appropriate scales. 




0.228 0.23 0.232 0.234 
sin 2 9 W (M z ) 



Figure 7: Scatter plot for {sin 2 9w(Mz), a 3 (Mz)}, as in Fig. |6], except that the 
intermediate matter thresholds are now also included in the analysis. 

Finally, we examine possible string models that contain additional (3,2)i/6 
representations. These representations are obtained, for example, in the flipped 
SU(5) x U{1) model |57j and in the SO(6) x £0(4) models |58|, and arise from 



the additional 16 and 16 representations. These may also arise in the standard-like 
string models. For example, in the model of Ref. a change in the sign of the 

GSO-projection phase C ^ 4 produces an additional (3,2)!/ 6 representation from 

the sector h^. However, a single (3,2)i/6 representation (or even several such repre- 
sentations) is not sufficient to accommodate the low-energy data, and indeed at least 
one (3, 1)1/3 representation is needed. With one pair of each, we find 

4 x 10 13 GeV < M 32 < 1 x 10 15 GeV , (5.53) 

with the corresponding constraints 

3 x 10 6 GeV < (M 31 ) uppcr < 1.5 x 10 15 GeV 
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5 x 10 2 GeV < (M 3 i) lower < 2.6 x 10 11 GeV . (5.54) 

In all of these cases, it is an important observation that such ranges exist in 
which string-scale gauge coupling unification is consistent with low-energy data. It 
is of course obvious that the presence of extra matter can in principle resolve the 
discrepancy between the GUT and string unification scales. What is highly non- 
trivial, however, is that precisely the required sorts of extra states naturally appear 
in some of the realistic free-fermion models we have examined, with the necessary 
non-standard hypercharges to do the job. Indeed, the crucial representations, as we 
have seen, are those for which the one-loop beta-function coefficients 62 and 63 are 
rather large, while b\ is small. It is for this reason that this particular string-predicted 
extra matter is able to modify the running of the strong and electroweak couplings 
without substantially affecting the £7(1) coupling. 

5.7 General conclusions 



It is apparent from the analysis presented in this section that string gauge coupling 
unification imposes a strong constraint on the allowed string models. Models that 
would otherwise provide a very attractive low-energy phenomenology can be ruled 
out on the basis that gauge coupling unification at the string scale cannot be in 
agreement with the low-energy data. The model of Ref. |2T| is an example of such a 
model. In this model, the extra color triplets and electroweak doublets that appear all 
have quantum numbers that fit into 5 and 5 representations of SU(5). Consequently, 
as we have shown, the extra states beyond the MSSM in these models cannot bring 
the string unification prediction into agreement with experiment. We find it very 
encouraging that some otherwise very appealing string models can be ruled out on 
this basis. 

Perhaps even more importantly, however, there exist realistic string models in 
which there naturally appears additional matter that does not fit into 5 or 5 rep- 
resentations of SU(5). As we have shown, for such string models the hypothesis of 
gauge coupling unification at the string scale can be in agreement with the low-energy 
data. The model of Refs. fT7|, E3J is one explicit example of such a string model which 



contains all the needed representations to achieve string gauge coupling unification, 
in just the right combinations and with just the correct hypercharges. 

We have seen, then, that intermediate matter thresholds seem to be the only pos- 
sible way in which string-scale gauge coupling unification can be achieved in realistic 
level-one string models. By level-one string models, we mean string models in which 
the non-Abelian gauge group content of the Standard Model is realized through a 
level-one Kac-Moody algebra. While our conclusions might be modified for models in 
which the gauge group is realized by Kac-Moody algebras at higher level, no realistic 
three-generation models of this sort have been constructed to date. Consequently, 
the need for (and appearance of) an additional matter spectrum beyond the MSSM 
seems to be a prediction of realistic level-one string models. 
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It is remarkable that string theory, which predicts an unexpectedly high unifi- 
cation scale M string , in many cases also simultaneously predicts precisely the extra 
exotic particles needed to reconcile this higher scale with low-energy data. As we 
have shown, all of the other possible effects are not likely to bridge the gap between 
the MSSM and string unification scales. This is a profound conclusion, and may have 
important experimental consequences. 



6 Why are string-induced threshold corrections so small? 

In this section we address the general question of the size of the string-induced 
threshold corrections. Our results in Sects. 4 and 5 for the realistic heterotic free- 



fermionic models, as well as similar previous calculations for orbifold models [19, [40 



and Type-II superstring models J4"I| , all consistently indicate that generic string- 
induced threshold corrections are relatively small. This result is a priori surprising, 
given the infinite towers of massive states which could in principle affect the running 
of the gauge couplings. Indeed, we have already remarked that the calculation of 
the heavy string threshold corrections Aq in (|3.8|) is similar to that of the one-loop 
cosmological constant A, 

A = I ^Z(r) , (6.1) 

and typical values of this cosmological-constant integral ( |6.1|) for non-supersymmetric 
string models are found |59j to be in the range A ~ (9(10 2 ). Taking this as a typical 
scale for such one-loop string amplitudes, the question then arises as to why a different 
particular amplitude, namely the one-loop heavy string threshold correction, should 
be so highly suppressed. 

6.1 General argument 

This suppression of the string-induced threshold correction is independent of any 
particular class of string model-construction, and is in fact a general property of 
string theory. We therefore seek a model-independent explanation of the underlying 
reason for this suppression. 

We begin by focusing on the behavior of the modified partition function Bg(t) 
which serves as the string integrand in the calculation of the threshold correction 
Aq in (|3.8| ), and in particular let us compare it with the ordinary partition function 
Z(t) which plays the same role in the calculation of the cosmological constant A. In 
analogy to ( |3.19| ), we can always expand the partition function Z in the form 

Z = 7-2- 1 £a mn g m <f (6.2) 



so that, in analogy with ( |3.21| ), we have 

A = } 



m,n 



^ a mn Imn • (6-3) 

m,n 
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There are therefore two potential sources of difference between the values of A and A: 
the integrals themselves are different due to the different powers of r 2 , i.e., 2^ ^ /(^ 
and the coefficients may also be different, i.e., b mn ^ a mn . Now, while it is true that 
\l£n\ i s usually less than this is not a major source of difference for two reasons: 

first, \I$ n \ is numerically smaller than |/^| by only a few percent, and second, the 
signs of these integrals tend to alternate with different values of m and n. Thus, the 
major difference between the values of A and A must involve the relative behavior of 
the coefficients {a mn } and {b mn }. 

In order to pin down the essential difference in the behavior of these coefficients, 
let us first recall their physical interpretations. In the ordinary partition function 
Z, each coefficient a mn represents simply the net number of states with spacetime 
squared masses (M|, M£) = (m,n), where "net" refers to the difference between 
the numbers of spacetime bosonic and fermionic states. By contrast, in the modified 
partition functions B, the coefficients h mn represent the charges of these states relative 
to the gauge group in question. Thus, states only contribute to b mn if they carry a 
non-zero gauge charge. 

Now, the general behavior of the ordinary state degeneracies a mn in string theory 
is well-known (see, for example, Ref. f43fl). In particular, there are two generic fea- 



tures which concern us. The first is the appearance of so-called "unphysical tachyons" 
in non-supersymmetric heterotic string theories — i.e., states which contribute to co- 
efficients a mn with m + n < ( "tachyonic" ) but m 7^ n ("unphysical"). More specifi- 



cally, it can be shown [|3| that any non-supersymmetric heterotic string theory which 
contains gravitons will have ao,-i 7^ 0. (Note that this does not imply that the cor- 
responding spacetime spectrum contains physical tachyons; indeed, the absence of 
physical tachyons requires only that a nn = for all n < 0.) By contrast, physically 
sensible string models will never have charged unphysical tachyons with energy con- 
figuration (m,n) = (0, —1). This occurs because any state with energy configuration 
(m,n) = (0,-1) must arise as the vacuum state in a Neveu-Schwarz sector, and 
since such a vacuum state is necessarily a gauge singlet, it cannot carry a non-zero 
gauge charge. We therefore find that although ao,-i 7^ for all non-supersymmetric 
heterotic string models, we must have 60,-1 = regardless of whether spacetime 
supersymmetry is present. 

This is a crucial distinction, because the potential contributions from the (m, n) = 
(0, —1) unphysical tachyons are typically larger than those from any other state. 
Indeed, the integral Iq-i i s typically an order of magnitude greater than any other 
integral I$ n . This feature alone, therefore, is responsible for a sizable reduction in 
the total threshold correction. 

Of equal importance, however, is the second generic difference between the coef- 
ficients a mn and b mn : their behavior as m, n — > 00. As is well-known, the number of 
physical states in string theory grows exponentially with energy, so that 

a nn ~ An~ B e CVE as n -> 00 (6.4) 
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where A, B,C are constants. This exponential growth in the number of states is 
the famous Hagedorn phenomenon which signals the existence of either a maximum 
(Hagedorn) temperature, or string phase transition. Indeed, as discussed in [Q, the 
existence of such exponential growth in the number of string states at high energy 
is directly related, through modular invariance, to the existence of physical and/or 
unphysical tachyons at the low energy of the string spectrum, and is hence unavoid- 
able. Of course, this rapid growth in the number of physical states does not lead 
to a divergent cosmological constant A, for there is an even stronger corresponding 
suppression for the integrals 1^: 

J n« ~ e ~°' n asn^oo (6.5) 

where C is a positive constant. Nevertheless, the fact that the a nn grow so quickly 
opens up a significant range of values of energy n for which the contributions of 
massive states of energy n to the cosmological constant are still sizable. Thus, the 
cosmological constant receives important contributions not only from the unphysical 
tachyonic states (and the physical massless states), but also from the first several 
massive states. 

By contrast, for the threshold corrections, this second source of contributions is 
often removed as well, for in the case of certain threshold corrections it can be shown 
that the coefficients b nn will exhibit growth which is at most polynomial rather than 
exponential f43|. This observation, which is ultimately related to the modular prop- 



erties of the modified partition function B(r), will be discussed below. In particular, 
we shall see that this suppressed polynomial rate of growth occurs for those func- 
tions B(t) which arise solely from gauge-charge insertions of the form QlQm with 
L 7^ M, and for which there are no contributions from charged unphysical tachyons. 
Indeed, a quick scan of the intermediate results listed in Sect. 4 for the appropriate 
tachyon-free models verifies that the values of A for those cases with L ^ M are 
further suppressed by a sizable amount relative to those with L = M. Thus, in these 
cases, not only are the contributions from the unphysical tachyonic states absent (as 
are the contributions from the massless states), but even the contributions from the 
massive states are extraordinarily suppressed. Therefore, in these cases, there are 
two features which combine to produce the unusual suppression of the L ^ M string- 
induced threshold corrections relative to the corresponding cosmological constant. 
Indeed, for certain free-fermionic embeddings of the gauge group (such as that of the 
flipped SU(5) models), these L ^ M insertions are the only ones which are relevant. 
This double-suppression mechanism is then directly responsible for the diminished 
size of total relative threshold corrections such as Asu(5) — Au(i)- 

Let us now briefly discuss the cases in which charged unphysical tachyons do 
appear. In these cases, although there will be charged unphysical tachyons making 
contributions to the threshold corrections A@, we still must have &o,-i — (for the 
general reasons indicated above). Thus, the unphysical (m, n) = (0,-1) tachyons 
which would have led to the largest contributions to the threshold corrections must 
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still be absent. For example, the charged unphysical tachyons that end up contribut- 
ing to the threshold corrections in the non-supersymmetric SU(3) x 577(2) x 577(1) 
string model discussed in Sect. 4.5 are all of a more harmless variety, with energies 
m and n never more negative than —1/4: 

t 2 - 1 B(t) = (4g 3/4 + 120g 7/4 + ...)q' 1/4 + non - tachyonic . (6.6) 

Indeed, these charged unphysical tachyonic states do not contribute nearly as much 
as the (0,-1) states, since their corresponding integrals are highly suppressed, with 
73/4-1/4 ~ —1.44 x 10~ 3 as compared to 7^-1 ~ — H-04. Furthermore, even though 
these milder charged unphysical tachyonic states may lead to exponentially growing 
values of b mn even for the L 7^ M basis insertions, modular transformations can be 
used to show that the resulting rate of exponential growth is also highly suppressed; 
this suppression is ultimately due to the fact that these tachyons are never as tachy- 
onic as the (0,-1) tachyons which would otherwise control the rate of growth [fi3|l . 
Thus, even in the cases that charged unphysical tachyons appear, the above argu- 
ments remain intact, for its two main features, namely that 

bo-i = and b mn <C a mn as m, n — > 00 , (6.7) 

remain unaltered. Thus the string-induced threshold corrections continue to be sup- 
pressed in spite of the appearance of such unphysical tachyons. 

Finally, we remark that this type of analysis can also incorporate the one case in 
which string-theoretic threshold corrections are known to be large |33|, [35|: namely, 
as the values of particular moduli (such as radii of compactification) are taken to 
infinity. Indeed, we shall see below that taking such a limit changes the modular 
properties of the modified partition function B{r) in so fundamental a way that 
the above suppression mechanisms ( |6.7|) no longer apply. Thus, even in this case, 
the sizes of the threshold corrections can still be understood as a consequence of the 
modular properties of the modified partition functions B(r) and the behavior of their 
coefficients b mn . 

6.2 The modular properties of the modified partition function 

We now turn to the modular properties of the modified partition functions B(t), 
and the resulting behavior of their coefficients b mn as m, n — > 00. 

As discussed above, if a given string model is devoid of charged physical or un- 
physical tachyons, then the corresponding modified partition function B(t) must have 
coefficients b mn which always vanish for all m + n < 0, and which often grow at most 
polynomially as m,n — > 00. However, this is behavior for the coefficients is certainly 
unusual, given the expectations based on analyzing ordinary partition functions Z(t). 
Indeed, it can be shown that any partition function Z(r) with a mn = for all 



m + n < must in fact have all coefficients vanishing, so that Z = identically; 
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otherwise, if these unphysical tachyonic coefficients are non-zero, then the partition 
function Z must have coefficients a mn which grow exponentially. It turns out that 
this is a general theorem which applies to all partition functions Z which correspond 
to string theories in D > 2 uncompactified spacetime dimensions (or equivalently, 
for all modular functions with modular weights k < 0). The question then arises as 
to how the modified partition functions Bq{t) manage to evade these constraints, 
and survive to exhibit non-zero coefficients b mn at higher levels despite having no 
tachyonic contributions in many cases. Moreover, we also wish to determine the un- 
derlying reason why this growth is fundamentally different depending on the form 
of the QlQm insertion, polynomial if L ^ M, and exponential otherwise. We also 
wish to understand why the limit of large moduli changes these results, and allows 
the corresponding threshold corrections to grow large. 

To answer these questions, we must first recall some facts about modular func- 
tions. In general, modular functions fi(r) (such as conformal- field-theoretic charac- 
ters Xi) transform covariantly under modular transformations, 

fi = {ct + df £ My /,(r) , (6.8) 

where the exponent k is called the modular weight and where the matrix is a 
mixing matrix which represents the particular modular transformation in the space of 
functions f\. One then builds a full holomorphic/anti-holomorphic modular invariant 
function (such as a partition function) by combining two such sets of characters fi 
and gi, each with modular weight k, in the form 



Z(r,r) = r 2 k £ Nti 9i ( T ) fj(T) (6.9) 



ij 



where r 2 = Im r and where Nij is a matrix chosen to satisfy M^NM = N (where M 
and M are respectively the modular transformation representation matrices in the 
spaces of functions Qi and fi). For example, for the classes of free-fermionic string 
models we have been examining, our fundamental modular functions are the rj and 
functions which appear in (|3.11| ) and fl3.12|) , and from ( |3.11|) we see that the total 



combined holomorphic and anti-holomorphic functions fi and Qi take the general 
schematic forms 

/ ~ rr 24 e 22 

g ~ 7T 12 e 10 . (6.10) 

Since the rj and G functions are modular functions with weights k = 1/2, we see that 
the total partition functions Z for these models all have the total modular weight 
k = —1. Note that this is consistent with the factor of r 2 _1 which appears in (|3.11|). 



Indeed, the general relation between the total modular weight k of the partition 
function Z and the spacetime dimension of the corresponding string theory is 

k = 1 - D/2 . (6.11) 
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For the modified partition functions B, however, this is no longer the case, for the 
effect of the charge insertions into the trace is to increase the modular weight from 
k to k + 2. We can see this easily as follows. Let us first define the more general G 
function of two variables z and r, 



6 



ZT = 



E 



(z+0)(n+a) Jn+a) 2 /2 



(6.12) 



so that our usual functions of a single variable r can be obtained from these more 
general functions by projecting to z = 0: 



6 



a 
IP 



-2-niaf) 



e 



ZT. 



2 = 



(6.13) 



Then, in terms of this generalized function, the singly primed function 0' which 
results from the insertion of a single charge operator Q into the trace can be written 
as the result of a single derivative with respect to z: 



0' 



T 



-2-Kia 



{2m dz 1/3 



a 



ZT 



2=0 



(6.14) 



Likewise, the doubly primed function 0", which results from the insertion of two 
charge-operator insertions, can be obtained from 0(z|t) via a second derivative with 
respect to z: 



0" 



a 
IP 



m dT 



2 $_ 

V 2m J dz 2 
a' 



ZT, 



2=0 



(3 



T 



(6.15) 



Now, the general Q(z\t) functions have modular transformations of the form 







z 


CLT + b 


CT + d 


CT + d 



(ct + d) 1/2 exp 



11TCZ 

CT + d 



Q{z\t) 



(6.16) 



where we have neglected overall r- and ^-independent phases and mixing matrices. 
From this result [in particular the exponent of the (ct + d) factor], we easily see that 
the modular weight of each 0-function is k = 1/2. However, using (|6.16|) and taking a 
derivative with respect to z before the projection to z = 0, we find that 0' transforms 
just like except with exponent k + 1 = 3/2. Likewise, a second derivative with 
respect to z further increases the exponent tofc + 2 = 5/2. Now, recall from ( |3.7| ) that 
the modified partition functions B contain a total of two (helicity) charge insertions 
for the right-moving anti-holomorphic sector, and two corresponding (gauge) charge 
insertions for the left-moving holomorphic sector. Thus, the total modular weight of 
the modified partition function B is not the negative value k = — 1 which we would 
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have expected for four dimensional theories, but rather k' — k + 2 = +1. Note that 
this is also consistent with the two extra factors of r 2 which were inserted along with 
the charge insertions into the trace in ( |3.7| ). 

It is for this reason that B can be non-zero even though it may contain no con- 
tributions from physical or unphysical tachyons. Indeed, the theorem mentioned 
above, which would have forced such functions to vanish, holds only for functions 
with negative modular weights. By contrast, as discussed in Ref. |43[], the coeffi- 



cients for tachyon-free functions with positive modular weights are expected to grow 
polynomially. 

Given that the coefficients h mn are expected to grow polynomially for tachyon-free 
theories, we still must explain why this growth is in fact exponential for the cases of 
gauge charge insertions QlQm with L = M, or for the limit of large moduli. Indeed, 
such exponential growth is more typical for partition functions of negative modular 
weight. 

It turns out that for the cases of insertions QlQm with L = M, this behavior 
is caused by a modular anomaly which prevents the modified partition functions 
B(t) from being truly modular- invariant. We can see how this anomaly arises as 
follows. We have already shown above in (|6.14j ) that a single charge insertion Ql is 



equivalent to differentiation with respect to z followed by projection to z = 0, so that 
if a certain trace fi(r) without any charge insertion transforms modular-covariantly 
as in ( p.8|) , then the same trace with a single charge insertion f[{r) will transform 
modular covariantly with increased modular weight: 

fl (^) = (cr + d)^ £ Mij /J(r) . (6.17) 

However, this is not the case for the traces f" (r) with two identical charge insertions. 
As shown in ( |6.15| ), the insertion of two identical charges Ql is tantamount to two 



z-derivatives or a single r-derivative, yet these derivatives are not covariant with 
respect to modular transformations. Indeed, it is easy to see that if fi(r) transforms 
as in (|6.8|) under r — > r' = (ar + b)/ (cr + d), then 

j j 

(6-18) 

While the first term is of the proper covariant form, the second term is not. Rather, 
the true "covariant derivative" on the space of modular functions of weight k is 
actually 

here the contribution from the second term cancels the anomaly caused by the first. 
Thus, double gauge charge insertions of the form QlQm with L = M destroy the 
modular covariance of the holomorphic or left-moving sector, transforming a modular 
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function of weight k — — 1 into not only a modular function of weight k + 2 = +1, but 
also an anomaly term consisting of a modular function of the original weight k = — 1 
divided by r 2 . 

It is this feature which ultimately explains why exponential growth of the coeffi- 
cients b mn can occur for these double gauge-charge insertions. Recall that regardless 
of the particular gauge-charge insertion for the left-moving sector, there is always a 
double charge insertion for the right-moving sector — this is the helicity insertion 

2 

Q H . Thus, for the L = M cases, the anomaly term from the right-moving sector can 
combine with the anomaly term from the left-moving sector to recreate a modular 
function with equal left- and right-moving modular weights k = — 1. Indeed, such 
a term will appear divided by two powers of r 2 , just as required for reproducing a 
proper modular function of weight k = — 1. Thus, for the L = M cases, a remnant of 
the original k = — 1 behavior survives despite the charge insertions. It is this which 
is ultimately responsible for the exponential growth rate for the coefficients b mn in 
these cases. 

A similar mechanism occurs in the cases of models for which a modulus (such as 
a radius of compactification) is taken to infinity. In these cases, we are effectively 
changing the dimensionality of the theory, so that once again the modular weight k 
of the modified partition function is decreased below zero.[| Indeed, the act of taking 
a radius of compactification to infinity amounts to replacing 

99 1 1 j , 

(6.20) 



in the modified partition functions B(t). This occurs because the B-functions, which 
had previously represented the lattice sums over discrete momentum- and winding- 
mode vacua, now become continuous integrals which can be evaluated, yielding mere 
extra inverse factors of r 2 in the decompactification limit. Thus, the effective modular 
weight of the remaining function B(t) is decreased in this limit of large moduli, and 
the above arguments then lead to larger values of threshold corrections. 

2 

Finally, we remark that since the helicity insertion Q H always takes the form of 
a double-charge insertion (or equivalently a single r-derivative, with its associated 
modular anomaly), the modified partition functions B(t) are never modular- invariant 
in any case. However, an improved derivation of the string-theoretic gauge coupling 
threshold corrections which is manifestly modular-invariant has recently appeared in 



Ref. [38]. In this analysis, manifestly modular-invariant results are achieved by explic- 
itly introducing an infrared regulator (provided by introducing a curved spacetime), 
and by taking into account the back-reaction from both gauge and gravitational inter- 
actions. Moreover, these modifications result in expressions for threshold corrections 
which yield their absolute sizes, rather than merely their relative differences. These 
expressions also incorporate the additional contributions (such as those due to dila- 
ton tadpoles) which must be included in the cases of string models without spacetime 

* We thank I. Antoniadis for discussions on this point. 
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supersymmetry. We hope to repeat our threshold analysis using these expressions in 
a future work. 

7 Conclusions 

In this paper we examined in detail the problem of gauge coupling unification 
in realistic heterotic string models. The class of models that we studied are all 
constructed in the free-fermionic formulation, and are among the most realistic string 
models constructed to date. We believe that this fact is not accidental, but may reflect 
deeper properties of string compactification that are at present unknown. Indeed, 
the free-fermionic models are constructed at a highly symmetric point in the string 
compactification space, and the Z 2 x Z 2 orbifold structure that is realized in these 
realistic free-fermionic models may be deeply connected to the existence of only three 
light generations in nature. 

Despite their many attractive properties, however, the realistic free-fermionic 
models (and string theory in general) predict that the gauge couplings unify at the 
string scale, which is of approximately 5 x 10 17 GeV. This implies that if the string 
spectrum below the string scale is assumed to be that of the MSSM, then the string- 
predicted values of the low-energy parameters sin 2 9 W (M Z ) and a 3 (M z ) will be in 
disagreement with the experimentally observed values. We explicitly illustrated this 
disagreement by evaluating the associated renormalization group equations including 
two-loop, Yukawa-coupling, and scheme- conversion corrections. 

The question then arises as to whether this disagreement with the low-energy 
observables can be used to rule out this class of realistic string models, or whether 
other effects may arise to alter this conclusion. To answer this question, therefore, it 
is necessary to examine all of the possible effects that can modify this result. 

One possibility is that the tree-level string predictions may be modified by the 
one-loop heavy string threshold corrections due to the infinite tower of heavy string 
modes. We developed a method for evaluating these threshold corrections in the 
free-fermionic string models, and were able not only to evaluate the string threshold 
corrections to any desired accuracy, but also to perform various non-trivial consis- 
tency checks in our analysis. We evaluated these string threshold corrections within 
a range of realistic free-fermionic models, and in general we found that the string 
threshold corrections are small and cannot explain the disagreement with the experi- 
mentally observed values. In fact, in all the cases that we studied, we found that the 
string threshold corrections tend to elevate the string unification scale by approxi- 
mately 20%, and consequently enhance the disagreement with experiment. Moreover, 
we found that the string threshold corrections are, in general, not significantly af- 
fected by the choice of the gauge group, the existence of spacetime supersymmetry, or 
the presence of charged unphysical tachyons. We were also able to provide a model- 
independent argument which explains why such threshold corrections are naturally 
suppressed in string theory, except at points with large moduli. Our argument re- 
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lied on only the modular properties of the modified partition functions B(t) which 
enter the calculation of the threshold corrections, and hence should have validity be- 
yond the class of free-fermionic string models studied here. Hence, we conclude that 
string threshold corrections cannot resolve the disagreement with the experimentally 
observed values in these models. 

Other possible corrections may arise due to stringy modifications of the U(1) Y 
hypercharge normalizations, light SUSY thresholds, additional gauge structure at 
intermediate energy scales, or additional matter thresholds at intermediate energy 
scales. We found that smaller values of the U(1) Y normalization k\ can result in 
agreement with experiment, but we argued on general grounds that such smaller 
values of k\ are not possible in self-consistent string models. To analyze the possible 
contributions from light SUSY thresholds, we parametrized the low-energy SUSY 
spectrum assuming universal soft SUSY-breaking terms, and found that light SUSY 
thresholds also cannot bring the string predictions into agreement with experiment. 
This conclusion is unaltered even if the assumption of universality is relaxed. We also 
examined the possibility of additional gauge structure at intermediate energy scales, 
and showed that within the context of these realistic free-fermionic string models, 
this also cannot resolve the discrepancy. 

Finally, we examined the effect of intermediate matter thresholds. Because we 
were able to rule out each of the other effects in the realistic string models, only the 
presence of additional non-MSSM matter could possibly reconcile string-scale unifi- 
cation with low-energy data. In the realistic free-fermionic string models, additional 
matter beyond the MSSM generically appears in the massless spectrum. This matter 
takes the form of color triplets and electroweak doublets, in vector-like representa- 
tions. While some of these extra matter representations have the weak hypercharge 
assignments that are common in grand-unified theories, some have weak hypercharge 
assignments that are unique to the string models and do not arise in regular GUT's. 
The ultimate mass scales of these extra states will be determined by the renormal- 
izable and non-renormalizable terms in the superpotential, but in general the mass 
scale of the additional states will not be at the string scale. Therefore, in general 
it is in unjustified to assume that the spectrum below the string scale is that of the 
MSSM. We showed that in some models, these appearance of these additional states 
does not resolve the discrepancy between string-scale unification and low-energy data; 
hence these models can be ruled out. More interestingly, however, we found that cer- 
tain other realistic string models provide just the right combinations of extra matter 
representations, with just the right gauge quantum numbers, to allow string-scale 
unification to be consistent with the low-energy data. Indeed, within these models, 
we found that a significant window exists for the additional mass scales so that the 
predicted values of the low-energy parameters a stTong (Mz) and sin 2 6 W(Mz) will be 
in agreement with the experimentally observed values. In some of these models (e.g., 
the model of Refs. |T7|, |2B|), these extra states have a uniquely stringy origin, and may 
have profound experimental implications. It is then imperative to examine whether 
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the mass scales of these additional thresholds can be derived from the string models, 
and take the desired values. Such work is currently in progress, and will be reported 
in future publications. 
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Appendix A 



Here we provide the explicit Vj and defining parameters for each of the models 
we consider in the text. 

As we discussed in Sect. 2, four-dimensional string models in the free-fermionic 
construction are defined by a set of boundary conditions for the 64 worldsheet 
Majorana-Weyl fermions, as well as a set of phases which ultimately describe how the 
generalized GSO projections are to be performed in each sector. In the notation of 



Ref. ITTTJ which we shall use below, these boundary-condition vectors are denoted Vj, 
i = 0, iV (where N + 1 is the model-dependent number of vectors which are neces- 
sary); likewise, the phases are defined through an (N + 1) x (N + 1) matrix denoted 
kij. However, most of the models we analyze originally appeared in the literature 
in the different notation of Ref. [10]. In this latter notation, boundary-condition 

b, 



vectors are denoted bj and the generalized GSO phases are denoted C 



b, 



For 



completeness, we now give below the explicit mapping between these two notations: 



C 



b, 



-2V, 



exp 



2ni (v^ + Vj + ^-V.-V. 



(A.1) 



where denotes the first component of the vector (i.e., the component cor- 
responding to the worldsheet fermion ^ carrying spacetime Lorentz indices). Fur- 
thermore, inner products such as Vj • V, are defined with opposite signs in the two 
notations, so that 

bj-b, = -4Vi-Vj . (A.2) 

Finally, however, we point out that an important issue is the adoption of a self- 
consistent scheme for handling the zero-modes of real Ramond fermions. In the 
notation of Ref. 1TTJ , an explicit convention is established and we have adopted this 



convention in our calculations; this is the origin, for example, of the crucial phase 



contribution which appears in ( 3.13 ). Without this phase contribution, our 
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expressions for the partition function Z(t) would not have been modular- invariant. 
However, in the literature (and, in particular, in the papers in which these models 
first appeared), other conventions have been implicitly adopted for handling the zero- 
modes. Thus, it is necessary to take these changes of convention into account when 
translating between the two notations. Fortunately, for the case of all of the present 
models, these extra changes of convention can be absorbed through changes in the 
A^-phases [beyond those kij values implied by ( |A.1|) ]. It is only with these additional 
effective phase changes that the correct particle spectrum for each model is produced. 
Thus, in terms of the exact (Vj, ky) notation and zero- mode conventions of 



Ref. [PH| , the models which we have analyzed in this paper can be listed as follows. 
First, as discussed in Sect. 2, they all share the same first boundary-condition vectors 
which comprise the so-called NAHE set. In the present V- notation appropriate for 
64 Majorana-Weyl fermions, the NAHE set takes the following form: 



right-movers 



left- movers 



V 

v 2 
v 3 
v 4 



11111111111111111111 

11100100100100100100 
11100100010010010010 
11010010100100001001 
11001001001001100100 



11111111111111111111111111111111111111111111 

00000000000000000000000000000000000000000000 
00001010101011111100000000001111110000000000 
10100000010111111010000000001111101000000000 
01010101000011111001000000001111100100000000 



where the entry '1' is shorthand for ' — 1/2'. Here the Vo through V4 vectors corre- 
spond to the vectors {1, S, bi, b 2 , b 3 } described previously. 

The individual models are then each derived from this underlying NAHE structure 
through the addition of extra boundary-condition vectors and the specific choice 
of GSO phases ky. These are listed below. 



A.l Flipped SU(5) model 

This is the model we considered in Sect, 
boundary-condition vectors 



4.1. This model has the additional 





right-movers 


left-movers 


v 5 
v 6 
v 7 


11100100010001001010 
11001010100100001010 
00000000000011000011 


00001001011011111100000000001111110000000000 
01100000011011111010000000001111101000000000 

000101110011 1100 1100 



where '— ' is shorthand for ' — 1/4'. Furthermore, fcoo an d the following phases k^ with 
i > j are equal to 1/2: k 2Q , k 2 x, k 30 , k 31 , k 32 , h , k 41 , k 42 , h 3 , k 51 , k 52 , k 53 , fc 54 , fc 61 , 
^62, k 63 , k 6 4, k 65 , k 70 , k 76 . The remaining phases with i > j are vanishing, and those 
with i < j are determined from those with i > j using the constraint equations given 
in Ref. 

Note that starting from this model as it was originally published and following 
the strict mapping given in (|A.l ), we would have instead found that ^50 and k^o 
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should also be non-zero (as well as corresponding changes in those i < j phases 



which are related to these through the constraint equations of Ref. [|ll|]). However, as 
discussed above, we must set k 50 and /c 60 to zero in order to account for the difference 
in the implicit conventions for dealing with the Ramond zero-modes. Only with such 
modifications is the correct model with the same particle spectrum reproduced. 



A. 2 First SU(3) x SU(2) x U(l) model 

This is the model we considered in Sect. 4.2. 
boundary-condition vectors 



This model has the additional 





right-movers 


left-movers 


v 5 
v 6 
v 7 


00000000011000000011 
00011000000000011000 
00000011000011000000 


10010101010111100000111100001110000011110000 
01011001101011100000111100001110000011110000 
101010100110 011—0 011—0 



where again '— ' is shorthand for ' — 1/4'. In this model, koo and the following phases 
kij with i > j are equal to 1/2: k 20 , k 2 i, k 30 , k 31 , k 32 , h , k 41 , k 42 , fc 43 , fc 50 , k 52 , k 60 , 
^63? ^65? ^74? ^75- All remaining i > j phases vanish, and those with i < j may be 



determined from those with i > j using the constraint equations given in Ref. ||11|| . 

In this case, no phase adjustments are necessary in order to account for the 
difference in Ramond zero-mode conventions. 



A. 3 Non-supersymmetric version of first SU(3) x SU(2) x U(l) model 

This is the model we considered in Sect. 4.5. This model has the same additional 
boundary-condition vectors and phases as the model above, except that fc 51 , fc 61 , 
and k 7 i are changed from to 1/2. (Of course, given the constraint equations of 
Ref. |ll], this induces corresponding changes in the coefficients k^ with i < j.) As 
required, these changes have the effect of breaking spacetime supersymmetry without 
introducing physical tachyons into the spectrum. 

A.4 Second SU(3) x SU(2) x U(l) model 



This is the model we considered in Sect. 4.3. This model has the additional 
boundary-condition vectors 





right-movers 


left-movers 


v 5 
v 6 
v 7 


00011011011011011011 
00000011000011000000 
00011000000000011000 


01100110100111100000111100001110000011110000 
10101010011011100000111100001110000011110000 
010110011010 011—0 011—0 



In this model, k 00 and the following phases k^ with i > j are equal to 1/2: k 2 o, k 2 i, 
&30, k 31 , k 32 , k 40 , k 41 , k 42 , k 43 , k 50 , k 52 , k 53 , fc 5 4, h , fc 65 , k 7Q , k 73 . All remaining i > j 
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phases vanish, and those with i < j may be determined from those with i > j using 



the constraint equations given in Ref. [11 



For this model, it is necessary to set k^ to zero in order be consistent with our 
Ramond zero-mode conventions. 



A.5 SO(6) x SO(4) model 

This is the model we considered in Sect. 4.4. 
boundary-condition vectors 



right-movers 



This model has the additional 



left- movers 



V 5 
V 6 

v 7 



11100100010001001010 
11001010100100001010 
00000000000000000011 
00000000000011000011 



00001001011011111100000000001111110000000000 
01100000011011111010000000001111101000000000 
00000000001111111111111100001111111111110000 
00010111001111100110000110001110011000011000 



In this model, fcoo and the following phases kij with i > j are equal to 1/2: /c 20 , k 2 \, 

^30) &31> ^32; ^40) ^41; &42; &43> &50 j ^51; ^52; ^53; ^60? ^65) ^70) ^72? ^73) ^74; 

fc 8 4, &86- All remaining i > j phases vanish, and those with i < j may be determined 



from those with i > j using the constraint equations given in Ref. [11 



For this model, in order be consistent with our Ramond zero-mode conventions, 
it is necessary to set /c 63 and /c 83 to zero, and k 53 and k$ 6 to 1/2. 
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